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The  thermostructural  response  of  cored  ceramic  bricks  used  in 
hypersonic  wind  tunnel  heaters  was  analyzed.  The  objective  was  to  gain 
an  understanding  of  the  mechanisms  which  can  cause  brick  fractures.  The 
analyses  were  carried  out  by  computing  (1)  the  gross  heater  thermal 
response  which  furnished  boundary  conditions  for  computation  of  (2)  the 
detailed  temperature  distribution  in  a  brick  seqment  which  was  input  to 
compute  (3)  the  stresses  and  strains  in  the  brick  segment.  State-of- 
the-art  nume-ical  analytical  techniques  were  applied  m  each  step. 

Yttria  stabilized  zirconia  and  alumina  bricks  were  considered.  Tt  was 
found  that  brick  "web  stresses"  (caused  by  temperature  gradients  normal 
to  the  hole  surface  during  heater  blowdown)  alone  are  not  severe  enenh 
to  cause  fractures.  Nonlinear  temperature  gradients  across  the  whole 
bricks  car  cause  significant  "body  stresses",  and  these  were  analyzed 
in  a  parametric  fashion.  The  effects  of  point  contacts  between  brick.' 
and  other  stress  affecting  phenomena  were  also  considered.  Super¬ 
position  was  applied  to  define  typical  combinations  of  thermal  and 
mechanical  loadings  which  could  cause  brick  stresses  equal  to  t-he 
estimated  brick  ultimate  strength. 
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ABSTRACT 

The  thermostructural  response  of  cored  ceramic  bricks  used  in  hypersonic 
wind  tunnel  heaters  was  analyzed.  The  objective  was  to  gain  an  understanding 
of  the  mechanisms  which  can  cause  brick  fractures.  The  analyses  were  carried 
out  by  computing  (1)  the  gross  heater  thermal  response  which  furnished  boundary 
conditions  for  computation  of  (2)  the  detailed  temperature  distribution  in  a 
brick  segment  which  was  input  to  compute  (3)  the  stresses  and  strains  in  the 
brick  segment.  State-of-the-art  numerical  analytical  techniques  were  applied 
in  each  step.  Yttria  stabilized'  zirconia  and  alumina  bricks  were  considered. 

It  was  found  that  brick  "web  stresses"  (caused  by  temperature  gradients  normal 
to  the  h^-le  surface  during  heater  blowdown)  alone  are  not  severe  enough  to  cause 
fractures.  Nonlinear  temperature  gradients  across  the  whole  bricks  can  cause 
significant  "oody  stresses",  and  these  were  analyzed  in  a  parametric  fashion. 

The  effects  of  point  contacts  between  bricks  and  other  stress  affecting  phenomena 
were  also  considered.  Superposition  was  applied  to  define  typical  combinations 
of  thermal  and  mechanical  loadings  which  could  cause  brick  stresses  equal  to  the 
estimated  brick  ultimate  strength. 
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SECTION  1 


INTRODUCTION 

Blowdown  type  hypersonic  wind  tunnels  utilize  storage  heaters  to 
increase  the  total  energy  content  of  the  air  used  as  the  test  fluid.  The 
maximum  operating  total  temperature  capability  is  usually  limited  by  the  high 
temperature  integrity  of  the  heater  matrix  material.  For  this  reason,  refractory 
ceramics  (e.g.,  alumina,  magnesia,  zirconia)  are  usually  employed. 

A  performance  requirement  for  hypersonic  tunnels  is  minimum  test  stream 
contamination.  The  primary  source  of  contamination  is  "dust"  resulting  from 
fractures  of  the  heater  core  matrix  material.  Thus,  the  requirement  to  maintain 
heater  matrix  thermal  stresses  safely  below  ultimate  strengths  results  in  major 
performance  and  design  constraints  on  the  wind  tunnel  facility  (e.g.,  Ref.  1). 

Most  existing  hypersonic  wind  tunnel  storage  heater  ceramic  matrices  are 
of  the  "pebble  bed"  type.  The  desire  to  eliminate  the  high  stresses  associated 
with  the  point  contact  of  stacked  spheres  has  led  to  the  identification  of  cored 
bricks  as  the  optimum  matrix  configuration  for  large  capacity  storage  heaters 
(Refs.  1,  2).  Cored  bricks  are  also  advantageous  in  that  they  enable  a  lower 
heater  pressure  drop  for  a  given  heat  transfer  area  density. 

Unfortunately,  the  t’..«.rraostructural  behavior  of  cored  bricks  in  wind 
tunnel  heaters  is  quite  complex  and  hence  poorly  understood.  However,  a  quanti¬ 
tative  understanding  of  the  thermostructual  response  of  the  heater  matrix  is 
necessary  to  enable  rational  heater  design  and  operating  envelope  definition. 
Cored  brick  thermal  stress  calculations  performed  to  date  have  been  quite 
approximate;  usually  involving  gross  assumptions  to  enable  the  application  of 
closed-form  solutions  or  handbook  formulas.  The  lack  of  understanding  of  the 
thermostructural  response  of  cored  bricks  is  evidenced  by  the  fracture  of 
bricks  used  in  pilot  heater  applications,  although  it  should  be  stated  that 
major  uncertainties  are  associated  with  factors  such  as  zirconia  destabilization 
(see,  e.g. ,  Ref.  3) . 

The  objective  of  this  study  was  to  apply  state-of-the-art  numerical 
analytical  techniques  to  analyze  the  thermostructural  state  of  cored  bricks  used 
in  hypersonic  wind  tunnel  heat  are.  The  analytical  techniques  included  multi¬ 
dimensional  numerical  transient  heat  transfer  analysis  codes  and  finite-element 
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structural  analysis  codes  previously  developed  primarily  for  application  to 
reentry  vehicle  he»*t  shields  and  nose  tips.  It  was  felt  that  understanding 
of  the  thermostructural  response  of  cored  brick  heater  matrices  could  be 
advanced  through  the  application  of  these  sophisticated  analytical  techniques. 

A  typical  configuration  was  selected  to  serve  as  a  basis  for  the  brick 
response  calculations.  Initially,  the  analyses  were  to  be  performed  for  conditions 
corresponding  to  a  pilot  heater  for  the  TTT  (True  Temperature  Tunnel)  facility 
once  planned  for  installation  at  AEDC.  However,  the  construction  of  this  facility 
remained  uncertain  and  so  conditions  corresponding  to  a  cored  brick  heater  planned 
for  the  AEDC  APTU  (Airbreathing  Propulsion  Test  Unit)  facility  were  selected  as 
a  basis  for  analysis. 

The  cored  bricks  are  of  hexagonal  cross-section  measuring  about  2.78  inches 
across  the  flats  as  sketched  in  Figure  1.  The  bricks  contain  91  holes  of  nominal 
0.2  inch  diameter.  This  results  i»  a  hole  spacing  to  diameter  ratio  (S/D)  of 
1.47.  Variations  on  this  nominal  hole  configuration  were  of  interest,  and  so 
S/D  =  1.3  and  1.8  were  also  considered.  The  heater  bed  is  86  inches  in  diameter 
by  24  feet  high.  The  matrix  consists  of  yttria  stabilized  zirconia  (YSZ)  bricks 
from  the  top  down  to  12  feet  (middle  of  the  bed)  and  alumina  bricks  from  the 
12  foot  level  down  to  the  bottom. 

The  nominal  operating  condition  considered  was  3000  psi,  1150  lbs/sec 
flow  rate,  with  a  bed  initial  temperature  distribution  of  sinusoidal  uhape 
given  by 

T(z)  -  Tm  -  1  (Tm  -  T0)jl  -  cos  [f.  (24  -  s) 

where  T  =  500 °R 

o 

T  «  4500nR 
m 

L  =  36  ft 

z  =  the  distance  from  the  bottom  of  the  heater  bed  in  feet 

Since  the  bed  is  24  feet  high,  the  initial  bed  temperature  is  45CC',R  at  the  top 
and  1500°”  at  the  bottom.  The  locations  in  the  matrix  of  primary  interest  are 
the  top  of  the  bed  (YSZ  bricks) ,  middle  of  the  bed  (YSZ  -  alumina  interface) ,  and 
the  bottom  of  the  bed  (alumina  bricks) . 
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During  tunnel  operation  ("blowdown"),  high  pressure  air  enters  the  heater 
from  the  bottom,  is  heated  as  it  passes  up  through  the  matrix  holes,  and  exhausts 
from  the  top  and  continues  on  to  the  test  section.  As  the  air  passes  through  the 
brick  holes,  heat  is  transferred  from  the  brick  to  the  air.  This  gives  rise  to 
sharp  temperature  gradients  in  the  brick  webs  between  the  holes,  the  hole  surface 
temperature  being  lower  than  the  web  center  temperature.  These  temperature 
gradients  cause  thermal  stresses  in  the  brick  webs. 

It  was  initially  suspected  that  these  "web  stresses"  were  the  most  critical 
stresses  in  the  brick,  and  so  major  emphasis  was  to  be  given  to  their  characteri¬ 
zation  for  a  variety  of  conditions.  However,  as  will  be  discussed  in  detail 
subsequently,  initial  calculations  demonstrated  that  the  web  stresses  alone  were 
very  much  lower  than  what  was  believed  to  be  the  material  strength.  Also, 
observation  of  fractured  bricks  from  a  pilot  heater  and  discussions  with  parties 
with  field  experience  with  cored  ceramic  bricks  indicated  that  brick  thermal 
stresses  arise  from  a  variety  of  factors  in  addition  to  web  stresses.  Rather 
than  limit  consideration  to  web  stresses,  it  was  decided  to  perform  calculations 
to  investigate  some  of  these  other  stress  effects. 

The  various  effects  which  may  give  rise  to  brick  stresses  may  be  analyzed 
with  various  degrees  of  precision.  Some  effects  are  not  particularly  amenable 
to  analysis  because  the  associated  boundary  or  initial  conditions  are  not  known 
'e.g.,  temperature  gradients  across  the  whole  brick  and  residual  stresses),  and 
so  these  were  analyzed  in  a  parametric  fashion. 

The  heater  and  brick  thermal  analyses  are  discussed,  and  some  typical 
results  are  presented,  in  Section  2.  These  analyses  provided  temperature 
distributions  which  were  input  to  the  stress  analyses.  The  brick  structural 
analyses  are  discussed  in  Section  3.  A  variety  of  potential  failure  modes  are 
discussed  and  quantitatively  characterized.  Conclusions  are  drawn  and  recommend¬ 
ations  are  made  in  Section  4 .  Certain  analytical  details  are  included  as  appen¬ 
dices.  Appendix  A  presents  some  mathematical  details  relative  to  the  heat  conduc¬ 
tion  solution  aspect  of  the  gross  heater  thermal  response  analysis.  Appendix  B 
documents  an  analysis  of  the  potential  effects  on  the  heater  thermal  response  of 
misalignments  of  the  holes  between  adjacent  bricks.  Appendix  C  discusses  some 
limitations  on  brick  stress  and  fracture  conclusions  due  to  combined  stress  effects 
and  fracture  criteria  uncertainties. 
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SECTION  2 


HEATER  AND  BRICK  THERMAL  RESPONSE 

An  accurate  detailed  thermal  history  or  a  given  region  within  the  heater 
core  must  be  supplied  as  input  to  any  thermostructural  analysis.  In  order  to 
assure  the  validity  of  this  thermal  analysis,  the  mathematics  employed  should 
model  as  well  as  possible  the  following  conditions  encountered  in  engineering 
applications: 

1.  Heater  beds  composed  of  multiple  materials  with  "real"  properties; 
i.e.,  materials  whose  conductivities  are  finite  and  whose  properties  are 
temperature  dependent 

2.  Variation  in  heat  transfer  coefficient  with  air  temperature  and  velocity 

3.  Initial  axial  temperature  distributions  of  practical  significance. 

Numerical  analysis  techniques  were  chosen  over  closed  form  solutions  for 
the  following  reasons: 

•  Closed  form  solutions  do  not  adequately  model  practical  applications 
because  they  require  the  assumption  of  a  single  matrix  material  with 
constant  thermal  properties  (with  infinite  conductivity  in  the  radial 
direction  and  zero  conductivity  in  the  axial  direction) ,  constant  heat 
transfer  coefficient,  and  uniform  initial  temperature  distribution 

•  The  basic  computation  schemes  and  in  some  cases  computer  codes 
already  existed 

•  Numerical  schemes  are  easy  to  interface  with  numerical  stress  calculations. 

The  basis  of  the  numerical  techniques  used  are  covered  in  Section  2.1,  the 
thermal  properties  of  the  materials  analyzed  are  given  in  Section  2.2,  and  typical 
results  of  the  thermal  analyses  are  presented  in  Section  2.3. 

2 . 1  ANALYTICAL  TECHNIQUES 

Structural  analysis  of  a  given  heater  region  requires  an  accurate  and 
highly  detailed  in-depth  thermal  history  of  that  region.  Numerical  thermal 
analysis  of  the  entire  heater  core  on  the  detail  scale  necessary  for  subsequent 
thermostructural  analysis  is  prohibitively  expensive  when  considered  in  light 
of  the  fact  that  detailed  thermal  histories  are  not  required  for  the  entire  bed 
but  only  at  selected  regions  in  the  bed  where  it  is  desired  to  calculate  the 


thermostructural  response.  Furthermore,  the  overall  thermal  response  of  the 
heater  may  be  adequately  calculated  without  determining  the  detailed  in-depth 
response  of  the  bed. 

In  view  of  these  considerations  it  was  decided  to  divide  the  thermal 
analysis  into  two  parts:  the  first  being  a  gross  heater  thermal  response  analysis 
to  determine  the  hole  surface  boundary  conditions  for  the  second  detailed  analysis 
of  a  given  heater  brick  segment. 

The  gross  thermal  analysis  is  presented  in  Section  2.1.1  and  the  detailed 
two-dimensional  transient  response  is  discussed  in  Section  2.1.2. 

2.1.1  Gross  Heater  Thermal  Response 

The  objective  of  the  gross  heater  thermal  response  calculations  is  to 
determine  the  transient  heat  transfer  conditions  (heat  transfer  coefficient 
and  gas  enthalpy)  existing  throughout  the  heater  bed  at  various  times  during 
the  blowdown.  A  modification  of  the  PANDA  (Ref.  4)  one-dimensional  conduction 
computer  code  was  utilized  for  these  calculations.  The  modified  computer  code 
is  referred  to  as  the  TACH  (Thermal  Analysis  of  Ceramic  Heaters)  code  and  the 
modified  in-depth  conduction  solution  is  discussed  in  detail  in  Appendix  A. 

To  ecoi  nically  accomplish  the  objectives  of  the  overall  heater  analysis 
the  simplifyir  j  assumptions  are  made  that  there  is  no  axial  heat  flow  in  either 
the  bed  material  or  the  flowing  air.  The  former  assumption  is  justified  on 
the  basis  of  the  relatively  lew  thermal  conductivity  of  almcst  all  ceramic 
materials  and  the  fact  that,  for  heater  geometries  of  interest,  the  radial 
temperature  gradients  are  an  order  of  magnitude  greater  than  the  axial  gradients. 
The  assumption  that  axial  conduction  within  the  air  is  negligibly  small  is  quite 
good  considering  the  relative  energy  transfer  due  to  bulk  fluid  movement. 

Figure  2  shows  an  overhead  view  of  the  holes  in  a  typical  heater  brick. 

Also  shown  are  the  lines  of  symmetry  between  the  holes;  no  heat  flows  across  these 
lines  of  symmetry  if  heat  transfer  conditions  are  constant  from  one  hole  to  the 
next.  This  latter  assumption  implies  that  flow  conditions  are  the  same  in  all 
holes  at  any  given  axial  location  and  that  losses  at  tne  edge  of  the  bed  are 
negligible. 

Notice  that  cymmetry  considerations  effectively  divide  the  bed  into  a 
series  of  nested,  adiabatic,  hexagonal  cylinders.  The  conduction  heat  transfer 


Significant  hole-to-hole  flow  variations  can  result  from  misalignment  t.£  the 
holes  in  adjacent  bricks,  and  this  can  give  rise  to  net  temperature  gradients 
across  the  bricks.  An  approximate  analysis  of  this  circumstance  was  carried  Tut 
and  is  presented  as  Appendix  B  to  this  report. 
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within  these  adiabatic  hexagons  is  Ivo-dimensional  in  nature  and  as  was  previ¬ 
ously  discussed,  need  not  be  completely  modeled  for  the  gross  heater  thermal 
analysis.  Therefore  in  order  to  reduce  the  problem  to  a  one-dimensional  con¬ 
duction  situation,  while  still  conserving  heat  storage  capacity,  an  effective 
adiabatic  circle  is  assumed  whose  area  is  equal  to  that  of  the  adiabatic  hexa¬ 
gon  (see  Figure  2) .  The  radius  of  this  circle  is  given  by 


R 


eff  " 


.525  S 


(1) 


This  approach  has  the  advantage  that  net  energy  is  still  conserved  whilo  reducing 
computation  time  by  at  least  two  orders  of  magnitude. 

Utilizing  the  lumped  parameter  approach  the  bed  is  divided,  for  numerical 
analysis,  into  a  number  of  isothermal  axial  stations  of  length  AZ.  Thus  the 
temperature  at  any  radial  location  is  assumed  constant  over  the  length  of  the 
axial  station  (AZ) .  Figure  3  depicts  the  essentials  to  be  considered  for  the 
defined  lumped  parameter  numerical  element. 

The  overall  solution  procedure  followed  in  predicting  the  transient 
thermal  response  of  a  numerical  element  is  as  follows: 

a.  A  suitably  small  time  step  is  determined  (as  described  in  Appendix  A), 
over  the  span  of  this  time  step  all  rates,  material  properties,  and 
boundary  conditions  are  assumed  constant. 

b.  A  system  energy  balance  is  performed  for  the  assumed  constant  rates 
and  boundary  conditions. 

c.  New  material  properties  and  boundary  conditions  are  determined  and 
the  first  step  is  then  repeated,  etc. 

The  system  energy  balance  is  discussed  in  Section  2. 1.1.1,  and  check 
solutions  are  shown  in  Section  2. 1.1. 2. 


2. 1.1.1  System  Energy  Balance 

Three  e^f.rgy  conservation  equations  may  be  written  for  the  numerical  ele¬ 
ment  of  Figure  3;  they  are 

1.  A  net  energy  balance  on  the  flowing  air 

2.  A  surface  energy  transfer  rate  balance  on  the  shell/air  interface 

3.  An  in-depth  energy  balance  on  the  cylindrical  shell. 
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The  net  rate  at  wnich  energy  is  entering  the  air  flowing  over  the  axial 
station  is 


Q  -  »<Hout  -  nln> 


(2) 


where 


m  =  air  mass  flow  rate 


H^n  =  bulk  enthalpy  of  incoming  air 


Hout  =  ^ulk  enthalpy  of  outgoing  air 

For  the  first  axial  station  the  enthalpy  of  the  incoming  air  (H^n)  is  a  constant 
equal  to  the  supply  air  enthalpy  and  for  all  subsequent  down  stream  stations 
Hin*  as  a  ^unct^-on  time,  is  known  from  the  transient  solution  at  the  previous 
station. 

The  enthalpy  of  air  as  a  function  of  temperature  is  input  tabularly  to 
account  for  temperature  variations  in  specific  heat  capacity. 

The  flux  of  energy  between  the  flowing  air  and  the  solid  heater  is  given 
by 

(3) 


q  =  h(Tw  -  Td) 


where 


h  =  the  local  heat  transfer  coefficient 
Tw  »  the  solid  wall  temperature 
Td  -  the  gas  "driving"  temperature 

For  purposes  of  calculation  the  gas  driving  temperature  (TD)  io  set  equal  to  the 
average  of  the  incoming  and  outgoing  bulk  gas  temperature 


TD  - 


Tin  +  Tout 


(4) 


The  local  heat  transfer  coefficient  (h)  is  calculated  by  the  method 
recommended  by  Colburn  as  presented  in  Reference  5: 

>  k 


h  =  .0231  (-^-)  Re'J  Pr1/3 


(5) 
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where 


kg  =  thermal  conductivity  of  the  fluid 
D  =  hole  diameter 

ReD  “  =  Reynolds  number  based  on  hole  diameter 

Pr  =  Prandtl  number  of  the  fluid 

Equation  5  represents  a  correlation  for  fully  turbulent  (i.e.  ReD  >  10,000) 
and  fully  developed  (i.e.Z  >  50  D)  pipe  flow.  The  Reynolds  numbers  of 
interest  here  are  typically  greater  than  25,000  and  Z  *  50  D  occurs  at 
approximately  10  inches  from  the  bed  inlet,  and  thus  both  conditions  are 
well  satisfied.  During  the  course  of  calculation  the  heat  transfer  coefficient 
is  constantly  re-calculated  (at  the  conditions  existing  during  the  previous 
time  step)  to  account  for  fluid  property  variation  with  temperature.  This 
latter  variation  can  amount  to  over  a  50%  change  in  heat  transfer  coefficient 
over  the  length  of  the  bed. 

The  in-depth  energy  balance  on  the  cylindrical  shell  of  heater  brick 
material  is  accomplished  by  a  fully  implicit  finite  difference  approach  which 
is  described  in  detail  in  Appendix  A.  Using  this  fully  implicit  method  the 
energy  flux  conducted  from  the  surface  is  given  by 

-  q  =  a  +  0TW  (6) 

where  a  and  6  are  defined  in  Appendix  A.  The  constants  c  and  B  result  from  a 
first  pass  Gauss  reduction  of  the  implicit  finite  difference  tri-diagonal 
conduction  matrix,  and  are  functions  of  the  material  properties  and  the  in-depth 
temperature  profile  of  the  previous  time  step. 

Combining  the  surface  energy  balance  with  the  results  of  the  in-depth 
energy  balance  (Equations  (3)  and  (6))  the  following  is  obtained: 

/-a  -  BTd\ 

q  *  hi  -  1  *  function  (Tn)  (7) 

\  h  +  6  / 

Note  that  the  net  heat  rate  to  the  flowing  air  must  be  equal  to  the  surface  heat 
flux  times  the  surface  area  or 
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q  A 


surf 


(8) 
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The  system  energy  balance  is  accomplished  iteratively  as  follows: 

1.  Guess  Tq 

2.  Calculate  q  from  the  combined  in-depth  and  surface  energy  balances, 
Equation  (7) 

•  . 

3.  Calculate  the  net  heat  rate  from  Q  =  q  Asurf 

4.  Calculate  HQut  from  the  net  energy  balance,  Equation  (2),  knowing  H^n 

5.  Look  up  T^n  and  Tout  from  H^n  and  HQUt  in  the  enthalpy  temperature  table. 

6.  Calculate  a  new  TD  from 


7.  If  the  new  Tq  is  within  an  allowable  error  of  the  old  TD  (say  +  5  deg 
Rankin)  go  to  the  next  time  step;  if  not  return  to  step  1.  with  the  new 
Tq  as  a  new  guess. 

Using  this  method  the  transient  thermal  response  of  each  axial  station  in  the 
heater  bed  is  calculated. 

2. 1.1. 2  TACH  Check  Solutions 

Following  its  development  phase  TACH  was  run  on  a  battery  of  test  cases 
to  determine: 

1.  Reliability  of  its  one-dimensional  conduction  solution 

2.  Agreement  of  numerical  solution  with  a  closed-form  solution  for  simple 
conditions 

3.  Optimal  nodal  spacing  and  time  step  size. 

The  one-dimensional  conduction  solution  of  TACH  was  checked  by  inputting  constant 
material  properties  and  running  for  one  axial  station  only.  The  output  was  then 
compared  with  an  analytical  solution  from  the  Schneider  charts  (Ref.  6) .  Good 
agreement  was  obtained  (less  than  1%  error). 

TACH  performance  on  the  total  bed  was  compared  with  Hausen's  analytical  sol¬ 
ution  (Ref.  7)  which  assumes  infinite  conduction  normal  to  the  flow,  constant 
material  and  gas  properties,  and  uniform  initial  axial  temperature  distribution. 

The  results  of  this  comparison  are  shown  in  Figures  4  and  5.  Fiqure  4  presents 
comparisons  of  non-dimensional  wall  temperature  ratios  (Tw/Twi)  versus  non-dimen¬ 
sional  distance  (?)  at  various  normalized  times  (n)  where 
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Tw  =  Twall  ~  Tentering  gas 
Twi  “  Twall  initial  “  'entering  gas 


Figure  5  shows  similar  comparisons  for  gas  temperatures  where 


t  *  T  —  T 
g  gas  entering  gas 

Again,  good  agreement  was  obtained. 

The  optimization  of  nodal  spacing  and  time  step  size  was  accomplished 
by  varying  either  one  by  a  factor  of  two  and  then  comparing  with  previous  results. 
This  optimization  resulted  in  the  selection  of  twelve  axial  nodes  over  the  length 
of  the  24  foot  bed,  five  radial  nodes  and  a  minimum  time  step  of  0.2  seconds. 

The  check  solutions  and  optimizations  lend  confidence  to  the  conclusion 
that  the  TACH  code  adequately  models  the  gross  transient  thermal  response  of 
ceramic  heater  beds  composed  of  multiple  materials  with  variable  properties  sub¬ 
jected  to  arbitrary  initial  temperature  distributions  and  accounting  for  variations 
in  heat  transfer  coefficient. 


2.1.2  Brick  Segment  Two-Dimensional  Transient  Response 


The  detailed  two-dimensional  transient  thermal  response  analysis  of  indi¬ 
vidual  brick  segments  was  performed  using  a  modified  version  of  the  ASTHMA  (Ref.  8) 
(Axi-symmetric  Transient  Heating  and  Material  Ablation  Program)  program.  The 
computation  method  used  is  an  alternating  direction  implicit-explicit  method  with 
an  internally  computed  time  step.  The  heated  wall  boundary  condition  is  specified 
by  a  tabular  list  of  the  time  dependent  heat  transfer  coefficient  and  gas  enthalpy 
at  the  bed  location  of  interest.  Both  of  these  quantities  were  calculated  by  the 
TACH  code. 


The  necessary  ASTHMA  modifications  involved  deletion  of  the  material  abla¬ 
tion  portion  of  the  code  as  well  as  conversion  from  axisymmetric  to  plane  two-ui- 
mensional.  The  conversion  from  axisymmetric  to  plane  two-dimensional  was  accomp¬ 
lished  by  noting  that,  for  a  constant  thickness  body,  the  axisymmetric  conduction 
solution  approaches  the  planar  solution  as  the  radial  dimension  becomes  large 
with  respect  to  the  body  thickness.  Therefore  a  two-dimensional  planar  conduction 
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solution  may  be  accomplished  with  an  axisymmetric  code  by  simply  treating  it  as 
a  pseudo  axisymmetric  problem  with  a  sufficiently  large  radius  of  curvature.  The 
necessary  radius  of  curvature  was  determined  from  conduction  theory  such  that 
computed  heat  fluxes  were  well  within  0.1%  of  their  planar  values. 

Figure  6  shows  the  adiabatic  hexagon,  previously  determined  from 
symmetry  considerations  (see  for  example  Figure  2)  together  with  additional  lines 
of  symmetry  that  exist  within  the  hexagon.  The  shaded  area  indicates  the  basic 
symmetry  element  of  the  core  heater  brick?  it  was  on  this  basic  symmetry  element 
that  all  detailed  two-dimensional  conduction  solutions  for  web  stress  analyses 
were  performed. 

The  thermal  grid  developed  to  model  the  symmetry  element  for  the  nominal 
cored  brick  geometry  is  shown  in  Figure  ?.  The  ASTHMA  code,  supplied  with  a 
thermal  grid  of  the  region  of  interest,  material  properties,  and  time  dependent 
boundary  conditions,  was  employed  to  compute  the  brick  segment  temperature  dis¬ 
tributions  for  any  desired  bed  locations  or  flow  conditions.  The  ASTHMA  code 
provides  punched  card  output  of  temperature  distribution  data  for  input  to  thermal 
stress  analyses. 

2.2  MATERIAL  THERMAL  PROPERTIES 

As  discussed  in  Section  1,  the  two  materials  of  interest  were  yttria  stab- 
alized  zirconia  and  polycrystalline  alumina.  The  thermal  properties  of  these 
materials  are  required  over  the  operating  temperature  range  of  the  heater  bed 
(i.e.,  1500°R  to  4500°R) . 

The  thermal  properties  used  and  the  sources  from  which  they  were  obtained 
are  covered  in  Sections  2.2.1  and  2.2.2  for  YSZ  and  alumina  respectively. 

2.2.1  Yttria  Stabilized  Zirconia  Thermal  Properties 

Zirconia  is  a  highly  refractory  material  with  an  operating  temperature  in 
excess  of  4500°R.  The  major  drawback  of  zirconia  is  that  it  experiences  a  des¬ 
tructive  crystalline  structure  change  in  the  temperature  range  of  from  2000  to 
2400°R.  To  counter  this  structure  change,  zirconia  is  stabilized  by  addition  of 
oxides  such  as  calcia,  magnesia,  or  yttria  (see,  e.g..  Refs.  1,  3). 

The  thermal  properties  of  yttria  stabilized  zirconia  (YSZ)  were  measured 
by  Southern  Research  Institute  and  obtained  by  Aerotherm  through  References  9 
and  10.  The  YSZ  thermal  properties  used  in  analysis  are  summarized  in  Table  1. 

2.2.2  Alumina  Thermal  Properties 

Alumina  is  a  commonly  used  refractory  material  with  a  maximum  service 
temperature  of  about  3600°R.  The  properties  used  were  for  the  more  stable  alpha 
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TABLE  1 


THERMAL  PROPERTIES  OF 
ZIRCONIA  USED  IN 


YTTRIA  STABILIZED 
COMPUTATIONS 


Temperature 

(°R) 

Conductivity 

k 

(Btu/ft-8ec°R) 
x  104 

Specific 
Heat  Capacity 
C 

(Btu/lb°R) 

500 

3.26 

.111 

1210 

3.26 

.130 

1960 

3.26 

.150 

2460 

3.31 

.161 

2960 

3.37 

.171 

3460 

3.47 

.178 

3960 

3.58 

.186 

4460 

3.94 

.195 

4700 

4.11 

.199 

Density  *  345 

lb/ft3 

TABLE  2 

THERMAL  PROPERTIES  OF  POLYCRYSTALLINE 
a  -  ALUMINA  USED  IN  CALCULATIONS 


Temperature 

<°R) 

Conductivity 

k 

(Btu/ft-sec°R) 
x  10“ 

Specific 

Heat  Capacity 

C 

(Btu/lb°R) 

500 

5.83 

.170 

1000 

2.86 

.255 

1500 

1.81 

.280 

2000 

1.22 

.310 

2500 

.92 

.320 

3000 

.83 

.330 

Density  =  248  lb/ft3 


alumina  in  a  polycrystalline,  100%  theoretically  dense  form.  Thermal  properties 
were  obtained  from  Reference  11  and  are  summarized  in  Table  2. 

2.3  TYPICAL  RESULTS  OP  THERMAL  ANALYSES 

Results  of  the  TACH  calculations  of  gross  heater  thermal  response  are  pre¬ 
sented  as  plots  of  wall  or  gas  temperature  versus  axial  length  at  various  times 
during  blowdown.  Figures  8  and  9  present  such  plots  for  the  nominal  bed  operating 
conditions  (i.e.,  mass  flow  rate  »  1150  lb/sec  and  sinusoidal  initial  temperature 
distribution  and  S/D  *  1.47).  Notice  the  effect  of  the  sinusoidal  initial  tempera¬ 
ture  distribution  on  reducing  strong  temperature  gradients  in  both  axial  length 
and  time. 

Typical  output  of  the  detailed  transient  response  analysis  (ASTHMA)  is 
shown  as  an  isotherm  plot  for  the  basic  brick  web  symmetry  element.  Figure  10 
presents  such  an  isotherm  plot  for  an  element  near  the  mid-bed  location  in  the 
zirconia  material  (Z  *  13  ft)  at  10  seconds  after  beginning  of  blowdown  for  the 
nominal  operating  conditions. 

Brick  wall  temperature  distributions,  such  as  shown  in  Figure  10,  were 
input  to  the  thermal  stress  analysis  which  is  described  in  the  next  section. 
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SECTION  3 

HEATER  BRICK  STRUCTURAL  P"SPGNSE 

A  variety  of  phenomena  can  potentially  give  rise  to  stresses  in  cor cc. 
ceramic  bricks  in  hypersonic  wind  tunnel  heaters.  Some  can  be  analyzed  quavtita- 
tively  with  high  confidence,  some  can  be  analyzed  parametrically,  while  others 
can  only  be  considered  qualitatively.  These  stress-affecting  phenomena  are  -".is- 
cussed  and  analysed  here  in  Section  3.3.  However,  prior  to  detailed  discuss.,  on 
of  the  analyses,  the  material  mechanical  properties  used  are  presented  (in  Section 

3.1)  and  the  numerical  stress  analysis  computer  code  used  is  discussed  (in  Section 

3.2) .  Conclusions  and  recommendations  based  on  these  analyses  are  presented  in 
Section  4. 

3.1  MATERIAL  MECHANICAL  PROPERTIES 

The  thermal-mechanical  properties  of  Ytfcria  Stabilized  Zirconia  (YSZ)  were 
determined  experimentally  by  Southern  Research  Institute  (Reference  12) .  Since 
the  material  is  isotropic  (mechanically) ,  the  relevant  properties  needed  for  a 
thermal  stress  analysis  are  (a)  E,  the  elastic  modulus,  (b)  v,  Poisson's  ratio, 

(c)  a,  the  coefficient  of  linear  thermal  expansion  and  (d)  Y,  the  yield  stress 
which  in  this  case  is  the  ultimate  stress.  All  of  these  properties  are  required 
for  the  temperature  range  0  -  4000°F.  The  shear  modulus  (G)  and  shear  yield 
stresses  (k)  were  determined  from  the  isotropic  formulas 


G 


E 

2(1  +  v) 


and  k 


Y 

7? 


The  material  behaves  in  an  elastic-brittle  manner,  and  consequently  the  yield 
stress  is  the  ultimate  stress.  The  experimental  results  were  interpreted  by 
both  AF ML  and  ARL  in  conjunction  with  Aerotherm  and  Weiler  Research  to  determine 
the  most  probable  valves  to  use  in  the  thermal  stress  analyses.  The  results  of 
this  interpretation  are  shown  in  Figures  11a  and  11b.  Poisson's  ratio  (v)  was 
assumed  constant  over  the  temperature  range.  Experimental  data  points  are 


shown  by  encircled  points  ir.  these  figures.  The  data  shown  above  34CG-F  are 


extrapolations  since  they  are  needed  for  the  analysis.  It  is  understood  that 


the  melting  temperature  of  YSZ  is  approximately  4500°F  and  so  the  valves  shown 


above  this  temperature  are  fictitious  and  only  used  for  the  purposes  of  the 


thermal  stress  analysis.  Additionally,  the  coefficient  of  linear  thermal 
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expansion,  a,  was  determined  from  the  thermal  strain  data  by  a  ■  efc^/(T  - 
T  )  where  TQ  was  the  reference  temperature  <70°F).  All  of  these  properties  are 
listed  in  Table  3  as  functions  of  temperature. 

TABLE  3 

YTTRIA  STABILIZED  ZIRCONIA  MATERIAL  PROPERTIES 


Temp 

°F 

E 

psi  x  10~6 

V 

a 

in/in  x  10+s 

Yield  Stress 
psi  x  10“3 

70. 

26.0 

.30 

.4390 

12.5 

2000. 

14.8 

.30 

.5606 

11.0 

3000. 

9.0 

.30 

.6229 

8.0 

3400. 

6.8 

.30 

.6360 

6.8 

4000. 

3.3 

.30 

.6514 

5.0 

5000. 

.01 

.30 

.6800 

2.0 

3.2  FINITE  ELEMENT  THERMAL  STRESS  ANALYSIS  CODE 

The  different  models  to  be  discussed  in  subsequent  sections  were  analysed 
using  the  DOASIS  thermal-stress  analysis  computer  code  (Ref.  13).  This  is  a 
general  2-dimensional  (plane  and/or  axisymmetric)  code  which  can  analyse  solid 
structures  composed  of  temperature  dependent  materials  and  loaded  by  thermal 
gradients,  boundary  pressures  and  body  forces.  For  the  present  cases,  only 
linear  thermal  elastic  stress  analyses  were  performed. 

To  use  this  code,  the  body  is  first  subdivided  into  a  large  number  of  small, 
discrete  pieces,  called  finite  elements.  Relationships  are  written  describing 
the  state  of  displacement  within  these  elements  based  upon  the  element's  tempera¬ 
ture,  material  properties,  etc.  These  relationships  (linear  algebraic  equations) 
are  then  connected  mathematically  such  that  displacement  continuity  exists  between 
the  elements.  The  resulting  set  of  equations  are  solved  and  yield  the  displace¬ 
ments  of  the  corners  (nodes)  of  the  elements,  from  which  the  strain  and  thus  stress 
within  each  element  may  be  calculated.  Once  the  resulting  stresses  are  known 
for  the  body  in  question,  iso-stress  contour  plots  can  be  made  which  show  lines 
of  constant  stress  existing  within  the  body.  For  each  of  the  different  models 
analysed,  a  finite  element  mesh  was  generated  and  the  appropriate  boundary  con¬ 
ditions  for  these  models  were  applied.  A  stress  analysis  was  performed  and  the 
resulti  were  plotted.  The  finite  element  discretation  (mesh)  used  will  be  shown 
along  with  the  resulting  iso-stress  contour  plots  for  each  of  these  models,  where 
appropriate.  For  further  details  concerning  structural  finite  element  analysis 
techniques,  see  References  13,  14,  or  15. 
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3.3 


BRICK  STRESS  ANALYSES 


The  objective  of  the  thermal  stress  analyses  performed  is  to  lend  insight 
into  the  question,  "what  type  of  thermal  stresses  developed  in  the  cored  bricks 
could  give  rise  to  thermal  stress  failures,  i.e.,  cracks  induced  by  thermal  stress?" 
To  accomplish  this  objective,  it  was  originally  planned  to  examine  primarily  the 
web  stresses  existing  at  the  hole  surfaces  due  to  the  rapid  brick  cooling  around 
the  holes  during  blowdown.  Subsequently,  additional  types  of  thermal  stress 
analyses  were  performed  to  help  lend  insight  into  possible  stress  producing  con¬ 
ditions  which  could  lead  eventually  to  fracture. 

Consider  the  cross  section  of  a  hexagonal  cored  YSZ  brick  as  shown  in 
Figure  1.  The  brick  measures  approximately  2.78  in.  across  the  flats.  This 
hexagonal  brick  contains  91  holes  in  a  hexagonal  pattern.  The  section  between 
the  holes  is  called  the  web.  During  a  blowdown  condition,  air  passes  through 
the  holes,  heat  is  rapidly  transferred  from  the  web  material  and  thereby  creates 
thermal  gradients  throughout  the  bricks.  If  the  geometry  of  the  bricks  was 
perfect,  the  connections  between  bricks  resting  on  top  of  one  another  perfect, 
flow  conditions  were  uniform,  and  no  heat  was  transferred  from  the  edge  of  the 
bricks  or  bed,  then  the  thermal  gradients  induced  in  the  bricks  due  to  this 
rapid  cooling  would  be  uniform  and  regular.  However,  since  these  conditions  are 
not  true  in  practice,  irregularities  do  occur  ar.d  therefore  ideal  conditions  are 
not  realized.  Consequently,  concern  arises  for  other  stress  producing  conditions 
which  will  occur  due  to  these  irregularities,  such  as  concentrated  loads  due  to 
two  irregular  surfaces  resting  on  top  of  one  another. 

The  major  effect  of  irregularly  shaped  bricks  is  that  nonlinear  thermal 
gradients  are  induced  across  the  bricks,  and  these  give  rise  to  thermal  stresses. 

The  reason  for  this  is  as  follows:  If  a  body  (with  constant  material  properties) 
is  heated  uniformly,  the  whole  body  expands  uniformly  and  no  stresses  are  induced, 
assuming  that  the  body  is  not  constrained  in  any  manner.  It  is  also  true  that 
if  a  linear  temperature  gradient  is  imposed  on  a  body,  the  body  will  deform 
linearly  and  again  no  stresses  are  induced.  This  follows  from  thermal  elas¬ 
ticity  theory  (e.g,,  Ref.  17).  Thus,  that  portion  of  a  temperature  gradient  which 
induces  thermal  stresses  into  a  body  is  the  nonlinear  part.  In  this  case,  po¬ 
tential  effects  of  nonlinear  temperature  distributions  across  the  heater  bricks 
were  investigated,  in  addition  to  the  previously  discussed  simple  web  stresses. 

To  summarize,  after  reviewing  the  cored  brick  heater  matrix  construction 
and  operation,  three  different  stress  conditions  were  analyzed  quantitatively 
relative  to  their  potential  for  causing  stresses  large  enough  to  cause  brick 
fracture.  First,  the  thermal  gradients  induced  into  the  web  sections  of  the 
bricks  during  blowdown  create  thermal  stresses,  as  previously  discussed,  since 
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these  thermal  gradients  are  not  linear  through  the  brick  web.  Secondly,  non6 
linear  thermal  gradients  can  be  induced  across  a  whole  brick  due  to  one  side 
being  cooled  faster  than  the  other.  This  would  arise  if  some  of  the  holes  in 
one  side  of  the  brick  didn't  match  up  correctly  and  thus  constrict  the  air  flow¬ 
ing  through  these  holes  (see  Appendix  B) .  Consequently,  both  radial  (across 
the  bricks)  and  axial  (along  the  bricks)  nonlinear  temperature  gradients  could 
be  induced  creating  thermal  stresses.  Thirdly,  if  the  end  surfaces  of  the 
bricks  are  not  perfectly  flat,  then  when  the  bed  is  assembled  by  stacking  bricks 
on  top  of  one  another,  these  non-flat  surfaces  will  transmit  the  weight  of  the 
bricks  at  points  (small  areas  of  contact)  rather  than  uniformly  across  the 
whole  cross  section.  This  will  give  rise  to  what  is  called  contact  stresses. 
Obviously,  the  compressive  components  of  these  contact  stresses  will  be  high, 
but  there  will  also  be  tensile  stresses  induced,  which  are  the  ones  of  impor¬ 
tance  here. 

The  next  three  subsections  will  present  the  analyses  of  these  three  dif¬ 
ferent  stress  conditions.  Following  these  discussion,  the  combined  effects  of 
these  three  will  be  discussed,  and  fracture  conditions  will  be  illustrated. 

Since  these  three  conditions  are  not  the  only  situations  which  can  cause  stresses 
to  occur  in  the  bricks,  other  effects,  which  do  not  easily  lend  themselves  to 
analytical  investigation,  will  be  discussed  qualitatively  in  subsequent  subsec¬ 
tions. 


3.3.1  Web  Stress  Analyses 


In  general,  an  accurate,  detailed  stress  analysis  of  any  one  cored  brick 
would  result  from  a  three-dimensional  model,  accounting  for  all  irregularities 
in  geometry,  etc.  However,  for  the  present  case,  the  thermal  stress  response 
of  the  cored  bricks  will  behave  fairly  uniformly  at  each  axial  station  along 
the  length  of  the  brick,  where  deviations  from  this  uniformity  will  occur  only 
on  the  ends  of  the  brick.  It  was  therefore  decided  that  a  generalized  two- 
dimensional  plane  strain  or  plane  stress  analysis  of  a  typical  cross-section 
would  yield  results  fairly  typical  of  what  any  detailed  three-dimensional 
stress  analysis  would  produce  for  stresses  in  thj  web. section  of  the  brick. 


The  term  generalized  plane  stress  (or  strain)  means  that  the  stress  (or 
strain)  normal  to  the  two-dimensional  plane  of  analysis  is  known  and  given  for 
purposes  of  analysis.  For  the  analyses  in  question,  three  different  types  of 
generalized  plane  analyses  were  performed.  First,  if  the  brick  is  sufficiently 
long,  a  cross-section  in  the  middle  of  the  brick  will  not  displace  axially  due 
to  the  material  above  and  below  it  preventing  such  movement.  This  constitutes 
a  state  of  plane  strain  in  which  the  axial  component  of  strain  is  zero.  Second, 
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if  the  brick  is  sufficiently  thin,  the  axial  components  of  stress  on  the  top 
and  bottom  surfaces  are  zero  and  therefore  axial  stress  developed  through  the 
thickness  of  this  section  cannot  be  too  diffeient  from  zero.  Therefore,  a  state 
of  plane  stress  exists  in  which  the  axial  component  of  stress  is  zero.  Third, 
if  one  considers  a  typical  section  existing  midway  in  the  whole  bed,  then  12 
feet  of  bricks  above  the  cross-section  of  interest  will  exert  a  dead  weight 
loading  on  the  cross-section  of  inte’est.  This  dead  weight  may  be  considered 
as  a  specified  axial  stress  equal  to  the  reacting  pressure  given  by 

p  =  density  x  height  =  345  x  12  ft  t  144  ^2—  =  28.9  psi 

ft3  ft2 

i.e.,  the  axial  component  of  stress  =  28.9  psi  and  is  constant  across  the  cross- 
section  of  the  bed.  Therefore,  a  state  of  generalized  plane  stress  exists  in 
which  the  normal  (axial)  component  of  stress  is  known. 

Now,  when  examining  the  geometry  of  the  cross-section  shown  in  Figure  1, 
an  additional  simplifying  assumption  may  be  made,  namely,  symmetry.  Due  to  the 
large  number  of  uniformly  placed  holes  (in  a  hexagonal  pattern)  one  would  ex¬ 
pect  a  symmetrical  temperature  response  to  exist  in  the  webs  as  discussed  in 
Section  2.1.2.  similarly,  a  symmetrical  thermal  stress  response  will  exist 
except  near  the  sides  of  the  cross-section  where  the  symmetry  conditions  do  not 
hold.  However,  in  the  center  of  the  brick,  symmetry  conditions  do  hold  and  may 
be  used.  Therefore,  a  typical  symmetry  web  section  was  chosen  for  the  above 
mentioned  three  types  of  plane  analyses.  This  typical  section  was  shown  in  Fig¬ 
ure  6.  The  actual  geometry  of  the  symmetrical  section  is  determined  by  the 
spacing  (S)  and  diameter  (D)  of  the  hole  pattern.  The  boundary  conditions  used 
in  these  plane  thermal  stress  analyses  are  shown  in  Figure  12,  where  the  sym¬ 
bolic  rollers  indicate  that  points  along  the  three  cut  boundaries  may  move  along 
these  cut  edges  but  not  normal  to  them.  The  fourth  side  which  is  the  hole  sur¬ 
face  is  considered  free. 

The  predicted  temperature  distributions  were  translated  from  the  thermal 
analysis  grids  to  the  stress  analysis  grids  following  the  extension  technique 
described  in  Reference  16,  Appendix  E.  This  extension  involves  an  extrapola¬ 
tion  of  the  temperature  field  from  interior  points  to  the  three  interior  bound¬ 
aries  formed  by  the  lines  of  symmetry.  Since  these  three  boundaries  were  all 
assumed  insulated  in  the  thermal  analysis,  this  extension  technique  was  simpli¬ 
fied  and  proved  to  be  highly  accurate  when  performed. 

Cue  to  the  three  boundary  conditions  shown  in  Figure  12,  any  uniform 
rise  in  temperature  applied  to  the  symmetrical  section  will  causa  a  uniform 
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expansion  of  this  material,  and  thereby  induce  compressive  stresses  which  in 
reality  would  not  be  present.  It  is  well  known  in  the  theory  of  thermal  stress 
analysis,  that  a  free  body  will  contain  zero  stresses  due  to  a  constant  or  lin¬ 
ear  varying  thermal  gradient  as  explained  earlier.  Considering  the  brick  as  a 
whole,  being  free,  then  a  uniform  rise  in  temperature  will  create  zero  stresses 
throughout  the  brick.  To  duplicate  this  condition  in  the  symmetrical  section 
web  stress  analysis,  the  weighted  mean  temperature  was  found  via 

/  T  dA 

Tmean  "  j  <ja 

and  used  as  the  stress  free  reference  temperature  in  the  analysis. 

The  nominal  S/D  ratio  chosen  for  the  web  stress  analyses  was  S/D  =  1.475 
with  a  hole  diameter  of  D  =  0.20  inch.  These  were  the  values  used  for  determin¬ 
ing  the  geometry  of  the  symmetrical  section  shown  in  Figures  6  and  12. 

The  in-plane  stress  analysis  for  the  symmetrical  section  described  here 
was  calculated  for  three  cases.  They  are:  Case  I >  az  =  -28.9  psi  which  corre¬ 
sponds  to  the  dead  weight  of  12  feet  of  YSZ  bricks.  Case  II,  =  0  which  is 

the  plane  stress  case  and  Case  III,  e  *  0  which  is  the  plane  strain  case.  The 

z 

time  at  which  stresses  were  to  be  analyzed  for  these  cases  was  selected  with 
the  objective  of  considering  the  maximum  thermal  stresses.  The  temperature  dis¬ 
tribution  predictions  were  interrogated  to  determine  the  time  at  which  maximum 

(Tweb  center  '  Thole  surface5  occurs*  This  corresponds  closely  to  the  time  of 
maximum  thermal  stresses  and  it  occurred  at  about  10  seconds  for  the  S/D  =  1.47 
nominal  blowdown  condition. 

The  stress  results  for  Case  I  are  shown  as  iso-stress  contour  plots  (i.e., 
lines  of  constant  stress)  in  Figures  13a  through  13f.  The  finite  element  dis- 
cretation  and  temperature  distribution  for  this  nominal  blowdown  case  at  time  = 
10.0  seconds  are  also  shown.  The  nominal  blowdown  case  is  for  p  =  3000  psi, 
mass  flow  rate  =  1150  lbs/sec  and  a  sinusoidal  initial  temperature  distribution 
T  =  4500°R  to  1520°R  with  a  brick  geometry  of  S/D  =  1.47,  as  described  earlier 
in  this  report.  The  axial  location  of  these  in-plane  stress  analyses  is  at  the 
YSZ  -  AI2O2  interface  which  is  at  the  middle  of  a  24-foot  column  of  bricks. 

When  examining  Figures  13  a  through  13f,  one  notices  that  the  maximum 

tensile  stresses  occur  at  the  hole  surface  and  are  almost  in  the  hoop  tension 

direction  around  this  hole,  with  a  value  of  0jjOOp  =  3000  psi  (c.f..  Figure  13f, 

amax5 '  Th^s  wa8  as  expected.  The  other  interesting  stress  plot  is  of  the  shear 

stress  t  shown  in  Figure  13e.  Here  we  see  that  a  large  amount  of  shear  distor- 
xy 

tion  occurs  near  the  hole  surface  and  along  the  longer  of  the  two  cut  symmetrical 
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sides  which  are  perpendicular  to  the  hole  surface  (along  the  top  side  of  the 
symmetrical  section  shown  in  Figure  13e).  This  shearing  is  caused  by  the  fact 
that  the  hexagonal  web  section,  when  subdivided  into  symmetrical  sections  as 
shown,  does  not  have  the  same  shape  as  sectors  of  an  equivalent  cylinder.  In  an 
equivalent  cylinder  sector  model,  the  shear  stresses  would  be  zero.  In  the 
analysis  presented  here,  they  are  not  zero  but  rather  40  percent  of  the  magni- 
tude  of  the  maximum  stresses.  This  is  significant  when  considering  the  overall 
effects  of  the  stresses  produced  within  this  symmetrical  section. 

The  plane  stress  (az  *  0)  case  (Case  II)  produced  results  which  were  al¬ 
most  identical  to  this  generalized  plane  stress  (a  =  28.9  psi)  case  (Case  I) 

z 

and  so  were  not  plotted.  The  reason  for  the  similarity  is  that  an  axial  stress 
of  oz  =  -28.9  will  interact  through  the  coupling  expressed  by  Poisson's  ratio 
(v  =  0.30)  and  only  produce  stresses  in  the  plane  of  approximately  8.5  psi  which 
is  negligible  as  compared  to  3000  psi. 

The  plane  strain  (e^  =  0)  case,  denoted  Case  III,  produced  significant 
tensile  stresses  in  the  axial  direction  of  the  bricks.  The  iso-stress  contour 
plots  (i.e.,  lines  of  constant  stress)  are  shown  for  this  case  in  Figures  14a 
through  14e.  The  finite  element  discretation  and  temperature  distribution  are 
the  same  as  for  the  other  cases  and  are  shown  in  Figures  13a  and  13b. 

Figures  14a  through  14e  show  that  the  maximum  tensile  stress  of  4250  and 
4000  psi  occur  for  the  hoop  tension  direction  (Figure  14e)and  the  axial  direc¬ 
tion  of  the  bed,  i.e.,  the  component  c  (Figure  14c),  respectively.  The  total 
state  of  stress  along  the  hole  surface  is  triaxial  since  the  componsnt  crx 
(approximately  normal  to  the  hole  surface)  is  also  tensile,  varying  between 
100  psi  to  1000  psi  as  shown  in  Figure  14a.  Therefore,  this  triaxial  state  of 

stress  a  :a  :a  varies  from  0.25  :  1.0  :  1.0  to  0.025  :  1.0  s  1.0  along  this 
a  y  z 

hole  surface,  the  peak  b-^ing  centered  near  the  longer  of  the  two  sides  which 
are  perpendicular  to  this  hole  surface  (the  top  side  shown  in  these  figures) . 

The  axial  stress  (cz)  developed  for  this  plane  strain  case  is  caused  by 

the  fact  that  the  colder  material  at  the  hole  surface  location  wants  to  contract 

in  the  axial  direction.  However,  the  hotter  material  in  the  center  of  the  web 

resists  this  contraction,  thereby  causing  this  hole  surface  region  to  go  into 

tension.  If  the  brick  is  long  enough  in  this  axial  direction,  then  the  center 

section  of  the  brick  will  be  in  a  state  of  plane  strain  (e  =  0) .  Shorter 

z 

length  bricks  will  deviate  from  this  pure  plane  strain  case  and  as  .he  length 
of  the  brick  goes  to  zero,  and  a  plane  stress  case  develops.  The  exact  magni¬ 
tude  of  plane  strain  (or  plane  stress)  conditions  for  short  length  bricks 
(approximately  2"  in  axial  length)  would  have  to  be  determined  by  a  three-dimen¬ 
sional  analysis. 
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FIGURE  14  RESULTS  OF  WEB  STRESS  ANAYSIS  "-EASE  IE. 
OONPITIONS  (PLANE  STRAIN,  6*  «  0)  NOMINAL  BLOW- 
QOWN  OONPITIONS,  TIME  —  10  SEC  ,  S/ P  =  1.47 
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The  results  of  the  maximum  stresses  for  all  three  cases  are  presented  in 
Table  4.  The  location  where  these  maximums  occur  are  indicated  in  Figure  13a. 
The  interesting  thing  about  the  results  presented  in  Table  4  is  that  the  plane 
strain  case  (Case  III)  produced  the  maximum  tensile  stresses,  not  only  in  the 
plane  but  also  normal  to  the  plane  (axial  direction  in  the  bed) .  Both  the  plane 
stress  (Case  II)  and  the  generalized  plane  stress  (Case  I)  solutions  produced  a 
biaxial  state  of  stress  (in  the  plane)  which  had  a  magnitude  of  only  approxi¬ 
mately  75  percent  of  the  plane  strain  case. 


TABLE  4 

MAXIMUM  STRESS  RESULTS  OF  IN-PLANE  ANALYSES 


Case 

max 

a 

X 

max 

0 

y 

max 

T 

xy 

max 

normal 

max 

0 

max 

max 

"max 

I  (oz  =  -28.192) 

750 

2764 

-1210 

-28.192 

3041 

1492 

II  (oz  =  0.0) 

747 

2751 

-1205 

0 

3028 

i486 

III  (e  =  0.0) 
z 

1071 

3952 

-1730 

4229 

4346 

213) 

The  web  stresses  analyzed  in  this  section  will  be  considered  with  respect 
to  other  stresses  and  the  material  strength  in  Section  3.3.4. 

3.3.2  Nonlinear  Temperature  Gradients  Across  Bricks 

Another  source  of  thermal  stresses  which  could  exist  in  the  cored  bricks 
is  due  to  non-uniform  heating  conditions  which  give  rise  to  non-uniform  tempera¬ 
ture  gradients  within  the  brick  as  discussed  in  Section  3.3.  Again,  to  repre¬ 
sent  this  behavior,  a  three-dimensional  model  would  be  the  most  accurate  course 
to  follow,  but  the  degree  and  type  of  nonlinear  temperature  gradients  is  not 
known.  To  keep  the  thermal  stress  analyses  within  economical  bounds,  and  in 
view  of  the  fact  that  the  actual  nonlinear  temperature  gradients  are  unknown, 
an  equivalent  cylinder  two-dimensional  model  was  used  with  various  hypothetical 
nonlinear  temperature  gradients.  The  limitations  of  this  equivalent  cylinder 
approximation  are  further  described  in  Appendix  C. 

This  equivalent  cylinder  model  was  developed  based  upon  the  hypothesis 
that  the  total  volume  of  the  equivalent  cylinder  equaled  the  total  volume  of 
the  hexagonal  brick.  Since  the  height  of  the  hexagonal  brick  and  cylinder  are 
the  same,  this  reduces  to  finding  the  radius  of  the  equivalent  cylinder  which 
will  make  the  total  cross-sectional  areas  the  same.  This  same  type  of  geometry 
problem  (i.e.,  find  an  equivalent  cylinder  for  a  given  hexagon)  was  solved  in 
Section  2.1.1  and  is  illustrated  in  Figure  2  where  the  equivalent  cylinder 
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cross-section  is  indicated.  This  equivalent  cylinder  for  the  entire  brick  has 
a  radius  given  by  a  =  W. 

The  boundary  conditions  of  the  two-dimensional  model  representing  the 
equivalent  cylinder  are  shown  in  Figure  15.  Notice  points  along  the  centerline 
are  free  to  move  axially  but  are  prevented  from  moving  radially  as  indicated  by 
the  rollers.  The  free  brick  is  fixed  against  rigid  body  motion  by  pinning  the 
point  on  the  centerline  at  the  bottom  surface.  The  remaining  three  sides  are 
free.  Two  basic  types  of  hypothetical  nonlinear  temperature  profiles  were  con¬ 
sidered  in  the  equivalent  cylinder  analyses.  These  were  axial  gradients  and 
radial  gradients  as  shown  in  Figure  16.  Tne  radial  gradients  are  specified 
by  two  parameters,  TBage  which  is  the  base  temperature  and  AT  which  is  the  de¬ 
gree  of  nonlinearity.  Assuming  a  quadratic  form  for  this  radial  temperature 
profile,  and  a  cylinder  radius  "a",  then  the  radial  temperature  (for  all  axial 
locations)  is  given  by 


T!r>  '  TBase  + 


(all  z) 


The  axial  gradients  are  specified  by  three  parameters,  TBottom  and  ttop  wh*cl1 
are  the  temperatures  at  the  bottom  and  top  of  the  equivalent  cylinder  respec¬ 
tively  and  AT  which  is  the  degree  of  nonlinearity.  Assuming  an  axial  height  of 
"h",  then  the  quadratic  form  for  axial  gradients  becomes 


T(z) 


+  K  (tt  -  V 


+  4 


(all  r) 


When  the  temperature  of  the  bottom  and  top  are  the  same,  the  second  term  in  the 
above  formula  disappears. 

The  complete  array  of  cases  considered  in  the  equivalent  cylinder  analyses 
for  both  radial  and  axial  gradients  are  listed  in  Table  5. 

The  reason  that  quadratic  nonlinear  temperature  gradients  are  considered 
is  that  for  a  free  body,  as  previously  discussed,  a  uniform  or  linear  tempera¬ 
ture  gradient  produces  zero  stresses  within  the  context  of  linear  thermal  stress 
analysis  (i.e.,  when  the  mechanical  properties  are  constant  and  not  functions  of 
temperature) .  Since  the  whole  temperature  range  considered  for  each  case  was  at 
maximum  only  at  3.3  percent  change,  the  temperature  dependent  properties  could 
be  considered  effectively  constant  in  this  range.  In  the  actual  analyses,  they 
were  considered  temperature  dependent,  but  still  the  degree  of  thermal  stresses 
produced  by  this  dependence  was  small  as  compared  to  the  thermal  stresses  pro¬ 
duced  by  the  nonlinear  portion  of  the  temperature  gradients. 
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The  stress  distributions  produced  by  nonlinear  temperature  distributions 
were  calculated  but  not  plotted.  The  finite  element  model  is  graphically  illus¬ 
trated  in  Figure  17.  The  results  of  these  stress  analyses  are  summarized  for 
the  maximum  tensile  values  existing  throughout  the  brick  in  Table  5.  There  are 
three  principal  stress  components  listed  in  Table  5.  They  are  or  (radial  stress), 
o0  (hoop  stress)  and  oz  (axial  stress) .  The  location  within  the  brick  where  the 
maximum  tensile  values  of  these  three  stresses  occur  is  shown  in  Figure  17  for 
radial  and  axial  temperature  distributions.  All  three  stresses  are  listed  for 
both  locations.  The  direction  of  these  three  principal  stresses  is  also  shown 
in  Figure  17  on  a  typical  small  element.  The  degree  of  nonlinear  temperature 
distribution,  AT,  is  listed  in  Table  5,  where  its  meaning  was  fully  explained 
earlier. 

When  examining  the  results  for  radial  nonlinear  temperature  distribu¬ 
tions,  one  notices  that  the  maximum  stresses  occur  on  the  outside  surface  of 
the  brick,  midway  between  the  top  and  bottom.  Here,  a  biaxial  stress  field 
exists  where  omax  =  o!1151.  Since  the  elastic  modulus  E,  coefficient  of  linear 

2  u 

thermal  expansion  a  and  the  AT  are  known,  then  the  principal  stresses  may  be 
expressed  by 


where  6^  is  a  numerical  factor  depending  upon  the  radial  temperature  distribu¬ 
tion.  Using  the  results  of  Table  5,  we  find  for  the  principal  stresses  er,  oQ, 
oz  for  the  temperatures  T  =  2000,  3000°F  the  following  6-factors: 

TABLE  6 

AT  (RADIAL)  6-FACTORS 


To  compare  these  numerical  results  for  a  6-inch  brick  with  a  closed  form 
analytical  solution,  we  select  a  solution  for  an  infinite  long  cylinder  having 
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TABLE  5 


STRESSES  ASSOCIATED  WITH 

NONLINEAR  RADIAL  AND  AXIAL  TEMPERATURE  GRADIENTS 


Type 

of 

AT 

3rick 

Length 

AT 

T 

Base  or 

Center  of  Ends  I 

Midway  Along  Outside 

Tbotton/Ttop 

°r 

°z 

°e 

°r 

cz 

°e 

20 

2000 

306.9 

-41.0 

306.9 

-62.9 

1172.3 

1229.0 

R 

A 

40 

2000 

610.3 

-81.7 

610. .3 

-125.8 

2343.5 

2458.5 

D 

60 

2000 

910.0 

-22.1 

910.0 

-188.7 

3513.4 

3688.7 

I 

A 

6" 

20 

3000 

196.1 

-26.2 

196.1 

-40.3 

751.3 

788.0 

L 

40 

3000 

387.9 

-51.8 

387.9 

-80.6 

1499.1 

1574.1 

60 

3000 

575.4 

-77.0 

575.4 

-120.7 

2243.3 

2358.3 

20 

2000 

162.3 

-7.2 

162.3 

-7.5 

211.7 

62.7 

4G 

2000 

323.4 

-14.4 

323.4 

-14.9 

421.6 

124.9 

60 

2000 

483.3 

-21.5 

483.3 

-22.2 

629.5 

186.6 

0 

200C/2100 

20.9 

.2 

18.5 

0 

-2.3 

-.8 

20 

2000/2100 

178.5 

-5.9 

177.4 

-7.3 

207.1 

61.4 

40 

2000/2100 

335.0 

-11.9 

335.2 

-14.7 

414.7 

122.9 

A 

C  M 

60 

2000/2100 

490.3 

-17.9 

491.8 

-21.9 

620.3 

183.8 

X 

I 

b 

20 

3000 

104.1 

-4.6 

104.1 

-4.8 

135.2 

40.1 

A 

L 

40 

3000 

207.3 

-9.2 

207.3 

-9.4 

267.6 

79.4 

60 

3000 

309.6 

-13.7 

309.6 

-14.0 

397.3 

118.0 

0 

3000/3100 

13.7 

.1 

12.4 

0 

-.8 

-.3 

20 

3000/3100 

115.2 

-3.8 

114.5 

-4.7 

130.9 

38.7 

40 

3000/3100 

215.6 

-7.7 

215.8 

-9.2 

259.8 

77.0 

60 

3000/3i00 

315.1 

-11.5 

316.0 

-13.7 

385.9 

114.5 

20 

2000 

356.9 

-15.2 

356.9 

-16.2 

447.2 

114.7 

4” 

40 

2000 

711.7 

-30.3 

711.7 

-32.2 

890.6 

228.5 

60 

2000 

1064.4 

-45.4 

1064.4 

-48.1 

1330.2 

341.3 

20 

2000 

1035.6 

-33.9 

1035.6 

-49.3 

816.4 

20.7 

2" 

40 

2000 

2068.1 

-67.9 

2068.1 

-98.1 

1627.4 

41.6 

60 

2000 

3097.4 

-101.8 

3097.4 

-146.5 

2432.5 

62.6 

* 

ft 

I 

* 
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FIGURE  17  LOCATIONS  OF  MAXIMUM  TENSILE  STEESS  IN 
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a  radial  temperature  distribution.  Jsing  this  together  with  the  previously  de¬ 
fined  temperature  distribution,  T(r),  we  find  that  at  r  •*  a,  the  maximum  stresses 
are  in  the  form  presented  above  with  the  numerical  factors  8  being 


er  =  o.o 
ee  -  o.s 

8Z  =  2v  =  0.6  (v  =  0.3) 

e 

These  analytical  results  are  very  close  to  the  B.  numerical  results  ex- 

*  s 

cept  that  since  only  a  6-inch  brick  was  used,  the  numerical  6X  is  slightly  low. 
This  comparison  was  made  midway  along  the  outside  location  where  o0  and  az  are 
maximum  (the  superscript  S  on  the  8  -  factors  indicates  this  location) . 

When  examining  the  axial  nonlinear  temperature  distributions,  the  follow¬ 
ing  table  of  8-factors  result. 


TABLE  7 

AT  (AXIAL)  6-FACTORS 


Center  of  Ends 

Midway  Along  Outside 

Length 

Temp 

•£ 

fiE 

6e 

4 

4 

8S 

e0 

2000 

.0682 

-.0030 

.0682 

-.0031 

.0889 

.0264 

2000/2100 

.0727 

-.0025 

.0726 

-.0031 

.0886 

.0263 

3000 

.0647 

-.0029 

.0647 

-.0029 

.0835 

.0248 

3000/3100 

.0702 

-.0025 

.0702 

-.0030 

.0833 

.0247 

4" 

2000 

.1501 

-.0064 

.1501 

-.0068 

.1878 

0.482 

2" 

2000 

.4362 

-.0143 

.4362 

-.0207 

.3432 

.0088 

No  simple  analytical  solution  was  found  in  the  literature  for  nonlinear 
axial  temperature  distributions  in  finite  length  cylinders,  and  consequently 
no  comparison  is  made.  One  interesting  fact  resulting  from  this  set  of  stress 
analyses  is  illustrated  in  Figure  18.  Here,  the  8-factors  for  both  the  radial 
and  axxai  nonlinear  temperature  distributions  are  shown  plotted  versus  the 
length  of  the  "equivalent  cylinder"  brick.  Acj  the  length  of  the  brick  decrease. ■ 
the  axial  8-factors  (and  therefore  the  stresses)  increase,  going  to  infinity  as 
the  length  goes  to  zero.  This  is  as  expected,  since  the  AT  was  held  constant 
(in  the  actual  stress  analyses)  as  the  length  was  shortened.  This  in  effect 
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concentrates  the  nonlinearity  into  a  short  length,  becoming  more  pronounced  as 
the  length  becomes  smaller  and  smaller  a  similar  effect  would  also  be  expected 
in  the  radial  0-factors  if  the  radius  was  decreased  while  holding  the  AT  constant. 

These  results  indicate  that  nonlinear  temperature  distributions  within  the 
bricks  can  produce  significant  tensile  stresses.  These  stresses  are  considered 
in  relation  to  other  stresses  and  brick  strength  in  Section  3.3.4 


3.3.3  Point  Contact  Loading  Stress  Analysis 


Another  type  of  loading  condition  which  could  give  rise  to  stresses  in 
the  bricks  is  non-thermal  in  nature.  This  is  concentrated  loading  conditions 
which  could  exist  due  to  the  bricks  transmitting  loads  from  one  to  another  at 
regions  of  small  area  instead  of  over  whole  load  transfer  areas.  One  such 
example  is  the  dead  weight  loading  of  bricks  stacked  on  top  of  one  another.  If 
the  ends  of  the  bricks  are  not  perfectly  flat,  then  this  dead  weight  loading  will 
be  transferred  at  points  rather  than  across  the  whole  end  surface  area.  This  is 
shown  graphically  in  Figure  19  where  one  end  of  a  brick  in  a  column  is  bowed, 
thereby  causing  the  dead  weight  of  all  the  bricks  in  the  column  to  be  transferred 
through  the  small  contact  area  which  here  is  the  outer  perimeter  of  the  brick. 


Again,  the  equivalent  cylinder  model  described  above  was  used  to  investi¬ 
gate  this  loading  condition.  The  model  and  boundary  conditions  were  the  same 
as  that  shown  in  Figure  15  with  the  exception  that  the  axial  constraint  at  the 
center  line  was  replaced  by  a  roller  condition  on  the  outside  bottom  circumfer¬ 
ence  which  allowed  radial  movement  but  no  axial  movement.  The  dead  weight 
loading  of  twelve  feet  of  bricks  (normal  pressure  of  28.9  psi)  was  applied  as  a 
normal  pressure  along  the  top  surface  of  the  equivalent  cylinder  model  of  the 
brick.  Force  equilibrium  then  produced  a  concentrated  line  load  along  the  out¬ 
side  bottom  circumference  where  the  axial  roller  boundary  condition  was  imposed, 
thereby  simulating  this  point  contact  load  transfer  condition. 


The  equivalent  cylinder  model  described  here  was  analyzed  for  this  mech¬ 
anical  line  loading  condition.  The  finite  element  model,  boundary  conditions 
(for  line-contact  loading)  and  loading  conditions  are  shown  in  Figure  20a.  The 
load  was  considered  as  the  dead  weight  of  a  twelve  foot  column  of  bricks  acting 
uniformly  over  the  top  of  the  six  inch  brick  model.  The  reacting  load  was 
considered  a  line-contact  load  existing  (axially)  on  the  outside  perimeter  as 
shown  by  the  support  condition  in  Figure  20a.  The  temperature  of  the  brick 
was  shown  by  the  support  condition  in  Figure  20a.  The  temperature  of  the  brick 
was  assumed  uniform  at  T  *  3000°F.  The  stress  results  are  shown  by  isostress 
contour  plots  in  Figures  20b  through  20f.  In  addition,  areas  of  tensile  and 
compressive  loading  conditions  are  indicated  on  these  figures. 
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These  stress  results  indicate  that  for  this  type  of  loading  condition, 

maximum  tensile  stresses  on  the  order  of  40  -  50  psi  are  developed  near  the 

area  of  the  line-contact  load.  This  may  be  seen  in  the  radial  stress  or  plot 

(Figure  20b) ,  hoop  stress  og  plot  (Figure  20d)  and  maximum  stress  plot 

(Figure  20f ) .  The  maximum  axial  stress  o  (Figure  20c)  developed  occurs  near 

z 

the  center  line  of  the  brick  and  is  approximately  2-3  psi.  In  contrast,  the 

maximum  compressive  o  is  approximately  150  psi  at  the  point  of  line-contact 

z 

loading.  The  maximum  shear  stress  Tmax  developed  is  approximately  95  psi  and 
also  occurs  at  the  point  of  contact  loading. 

In  an  overall  sense,  this  line-contact  loading  condition  does  not  produce 
any  significant  tensile  stresses  (only  40  -  50  psi)  as  compared  to  the  in-plane 
thermal  stresses  or  nonlinear  temperature  distribution  thermal  stresses.  How¬ 
ever,  the  stress  patterns  and  magnitudes  which  it  does  produce  can  contribute 
to  a  critical  condition  in  conjunction  with  the  other  modes  of  induced  stresses 
within  the  bricks. 

3.3.4  Combined  Effect  of  Different  Loadings 

When  examining  the  stresses  produced  by  the  three  different  sources  of 
loading,  none  of  these  loading  conditions  produces  sufficient  high  tensile 
stresses  to  cause  fracture  by  itself.  Since  the  material  acts  in  an  elastic- 
brittle  manner,  fracture  is  defined  here  as  a  tensile  stress  which  equals  or 
exceeds  the  yield  (or  ultimate)  allowable  stress  shown  in  Figure  11a.  When 
comparing  predicted  stresses  to  this  ultimate  stress,  the  combined  effect  of 
the  various  loadings  discussed  in  previous  sections  should  be  considered.  This 
can  be  done  through  superposition. 

As  an  illustration  of  superposing  two  thermal  stress  loadings,  consider 
the  stresses  produced  by  a  nonlinear  axial  and  radial  temperature  gradient. 

For  the  range  of  the  AT's  considered,  the  material  properties  may  be  considered 
effectively  constant.  Consequently,  the  principle  of  superposition  may  be  used 
to  combine  a  nonlinear  radial  temperature  stress  solution  and  »  nonlinear  axial 
temperature  stress  solution.  This  occurs  since  the  governing  equations  in 
thermal  elasticity  theory  are  linear  (when  the  material  properties  are  constant) . 
Now,  the  maximum  axial  stresses  in  the  equivalent  cylinder  brick  model  shown  in 
Figure  17  occurred  on  the  outside  surface  midway  along  the  brick  for  both  a 
radial  AT  and  an  axial  AT.  for  the  radial  at  loading,  the  maximum  hoop  stzess 
also  occurred  at  this  location  but  not  for  the  axial  AT  loading.  The  maximum 
hoop  stress  for  the  axial  AT  loa'  ing  occurred  at  the  top  (or  bottom)  of  the  model 
at  the  axis  of  revolution  of  the  model  (Figure  17).  However,  tensile  hoop 
stressed  do  exist  at  this  midway,  outside  location  and  are  given  in  Table  7  by 
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S  E 

the  -  factor  (as  compared  to  the  80  -  factor) .  Interpolating  from  Table  7 
for  a  temperature  of  T  r  2872°F  (the  stress  free  reference  temperature  used  in 
the  in-plane  solutions),  we  have  for  a  six  inch  brick: 

6?  =  .0842 
z 

(Axial  AT) 

8?  =  .0250 

V/ 

c 

When  comparing  the  8  -  factors  in  Table  7  for  six  inch  bricks  and  four  inch 
bricks,  one  notices  that  the  four  inch  brick  8  -  factors  are  approximately 
2.15  times  those  for  six  inch  bricks.  Hence,  assuming  this  generally  holds 

g 

true,  tne  above  8  -factors  for  a  four  inch  brick  becomes: 

6S  =  0.538 
z 

(Axial  AT) 

8g  =  .1810 

c 

Interpolating  from  Table  6  for  the  radial  AT  8  -  factors  at  the  temperature  T  = 
2872°F  yields: 


B?  =  .4900 
z 

B*  =  .4760 


(Radial  AT) 


Now,  the  total  axial  stress  cz  and  hoop  stress  at  this  outside,  midway 
location  due  to  both  radial  and  axial  / 7's  are  given  by: 

‘  T. &T  [-1810  “Axial  4  •4,6°  4Tn.diai] 

*  TT^T  f-0538  ‘-Axial  +  •490°  “Radial] 

Interpolating  for  the  material  properties  E,  a  and  v  from  Figure  11  for  a 
temperature  of  T  r  2872°F  yields: 

E  =  9 . 74  2x10 ‘psi 
a  =  6.149xiO"‘in/in°F 
v  =  0.300 
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we 


and  therefore: 


Eo 


=  85.88  psi/°F 
The  formulas  for  the  axial  and  hoop  stresses  now  become: 


°e  ’  4'604  “axial  4  40-74  “Radial 


c2  -  15.49  *  41.93 

These  formulas  represent  the  stresses  induced  in  a  four  inch  brick  at  T  * 

2872°F  at  the  outside  surface,  midway  in  the  axial  direction  of  the  brick, 

for  a  given  £TAx^ai  and  ^TRa(jiai*  For  example,  an  axial  AT  *=  148.  °F  and  a 

radial  AT  =  75. °F  would  yield  a  hoop  stress  of  cQ  =  3737.  psi  and  an  axial 

stress  of  a  =  5438.  psi. 

2 

Since  the  Zirconia  is  considered  an  elastic-brittle  material,  the 
Coulotnb-Mohr  failure  criteria  should  be  used  to  predict  fracture  conditions. 
In  the  present  case  where  oQ  and  oz  are  both  tensile  (neglecting  the  radial 
component  oR  and  the  shear  stress) ,  this  means  that  the  maximum  stress  a* 
due  to  oQ  and  oz  must  reach  the  ultimate  tensile  stress  condition  otJlt  to 
fracture,  i.e.. 


o*  *  =t>  Fracture. 

In  a  two  dimensional  stress  plane  (tensile  quadrant) ,  this  fracture  surface 
is  represented  by  a  circle  as  shown  in  Figure  21  (cf.  Reference  18).  For  the 
two  example  tensile  stresses  cQ  and  a z  given  here,  we  have: 

o*  =yfoJ  +  ai  =  6598.  psi 

v  Z 

The  ultimate  allowable  stress  at  T  =  2872 °F  is  given  by  the  uniaxial  data 
presented  in  Figure  11a  as  =  8384.  psi.  Therefore,  the  maximum  stress 

a*  =  6598.  psi  produced  by  a  ^TAxiai  =  148.  °F  and  a  ATRad^a^  “  ?5.*F  would  be 
78.7%  towards  fracture. 

This  illustration  of  superposition  only  considered  two  sources  of  loading, 
namely  a  radial  AT  and  an  axial  AT.  To  be  complete,  one  should  consider  all 
sources  of  loading  to  determine  maximum  stress  conditions.  Therefore,  the 
following  sources  should  be  considered : 
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I.  Body  thermal  stresses  due  to  nonlinear  axial  temperature  gradients. 

II.  Body  thermal  stresses  due  to  nonlinear  radial  temperature  gradients. 

III.  Web  thermal  stresses  (in-plane  solutions). 

IV.  Body  stresses  due  to  hydrostatic  pressure  of  p  =  3000  psi  acting  on  the 
whole  bed. 

V.  Local  stresses  due  to  point  or  line  contact  loading  conditions. 

VI.  Other  sources  of  stresses  such  as  local  hot  spots,  residual  stresses 
from  fabrication,  etc.  which  will  be  described  subsequently. 

The  last  two  sources  of  stresses  will  be  ignored  for  the  present  and 
only  the  first  four  will  be  considered. 

For  the  example  described  above  (i.e.,  due  to  sources  I  and  II),  the 
tensile  thermal  stresses  due  to  III  and  the  compressive  stresses  due  to  IV 
should  also  be  considered.  From  the  in-plane  web  thermal  stress  solutions 
described  earlier,  the  most  severe  case  of  web  thermal  stresses  due  to  blow¬ 
down  was  the  plane  strain  (cz  =  0)  solution,  i.e.,  Case  III.  From  Table  4, 
the  maximum  axial  and  hoop  stresses  of  Case  III  existed  at  the  same  location 
and  were: 


(c 


max 


5 hole  =  4346  psi 


(02>hole  =  4229  psi 


where  (a  ).  .  was  the  maximum  hoop  stress  existing  at  the  hole  surface  and 
ITldX  noi6 

(a  ).  ,  was  the  "normal"  or  axial  stress  from  this  in-plane  solution.  The 
z  hole 

(o  ),  ,  stresses  exist  in  the  plane  and  so  would  be  additive  to  the  o.  stresses 
max  hole  r  o 

of  the  equivalent  cylinder  solution.  The  (o  )j^e  stress  is  in  the  axial  direc¬ 
tion  of  the  brick  (normal  to  the  plane)  and  therefore  additive  to  the  o  stress 
of  the  equivalent  cylinder  solution.  The  stresses  due  to  IV  are  compressive, 
being  the  same  in  all  directions,  and  therefore: 


0^  =  -3000  psi 

a  ■*  -3000  psi 
z 

Now,  when  the  four  sources  of  stress,  i.e.,  I,  II,  III  and  IV  listed  above  are 
combined,  the  resultant  tensile  stresses  become: 
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o0  =  3737.  +  4346.  -  3000.  =  5083.  psi 


o  =  5438, 
2 


+  4229.  -  3000.  =  6667.  psi 


and  the  resultant  maximum  stress  o*  is: 

o*  =  Voo  +  o2  =  8384  .  psi 

which  is  just  the  fracture  stress  for  T  =  2872. °F.  Therefore,  when  con¬ 

sidering  the  first  four  sources  listed  above,  with  a  AT^^a^  =  148.06  and  a 
ATRadiai  =  75.0,  the  resultant  maximum  tensile  stress  is  just  sufficient  to 
fracture  the  brick  on  the  outside  surface,  midway  along  the  axial  direction 
of  the  brick.  It  is  also  possible  to  have  other  combinations  of  ATAx^a^  and 
ATRadial  combinin9  together  with  the  stresses  from  sources  III  and  IV  to  yield 
maximum  stress  which  equals  the  ultimate  stress  =  8384.  psi  for  this 

temperature  of  T  =  2872. °F 

The  formulas  for  the  oe  and  oz  are  given  by 

°e  =  4-604  “Axial  +  40-74  “Radial  +  4346'  ‘  3°00- 

“  15'49  “Axial  +  41'93  “Radial  4  4229'  '  30C°- 
and  the  fracture  condition  is  given  by 


°e  +  °z 


Juit  =  8384,  Psi’ 


These  three  equations  define  an  interaction  curve  between  AT  and  ATRad^ajL 

which  will  yield  a  fracture  condition.  This  interaction  curve  is  shown  plotted 
in  Figure  22.  The  values  of  ^TAx^a^  and  ^TRa(j^a^  used  in  the  example  given 
previously  are  shown  by  dashed  lines  in  Figure  22. 

The  reason  that  the  stresses  due  to  source  V  were  not  included  is  that 
the  contact  line  loading  solution  indicated  that  the  maximum  tensile  stresses 
developed  were  only  40  -  50  psi,  and  were  in  the  region  near  the  end  of  the  brick 
instead  of  near  the  middle.  However,  this  does  not  mean  that  a  point  contact 
loading  condition  near  the  middle  of  the  brick  (axially  speaking)  would  not 
contribute,  Lut  high  compressive  stresses  m  the  direction  of  the  point  (or  line) 
contact  load  accompany  the  relatively  small  tensile  stresses  developed  and  so 
this  source  of  tensile  stresses  was  neglected  as  a  major  contributor  to  a  condi¬ 
tion  which  would  produce  a  fracture. 

The  stresses  listed  in  source  VI  will  be  discussed  in  the  next  section. 
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3.3.5  Other  Stress-Affecting  Phenomena 

The  analyses  described  in  Sections  3.3.1  through  3.3.4  considered  only 
deterministic  phenomena  which  lent  themselves  to  analytical  investigation. 

Other  stress-affecting  phenomena  exist  within  the  bricks  which  cannot  be  treated 
analytically,  at  least  in  an  economical  manner.  These  other  phenomena  are  stress 
conditions  arising  from  (a)  built  in  residual  stresses,  (b)  local  hot  spots,  (c) 
partial  destabilization  of  the  zirconia,  (d)  micro-cracks  existing  in  bricks  and 
other  "real  world"  effects  which  are  usually  not  considered  in  analytical  pre¬ 
dictions.  Gome  of  these  phenomena  are  always  present  in  structural  components 
but  are  often  not  important  since  the  materials  are  ductile,  and  allow  local 
stress  rising  conditions  to  be  relieved.  However,  in  ceramic  materials  like 
yittria  stabilized  zirconia  (YSZ) ,  the  material  is  elastic  brittle,  and  does 
not  possess  any  ductile  behavior.  Stress  rising  conditions  will  not  be  relieved, 
and  will  be  the  mechanisms  which  initiate  fracture.  Tn  this  section,  some  of 
these  other  stress  rising  phenomena  are  discussed  in  a  qualitative  fashion. 

Residual  stresses  occur  during  the  ceramic  brick  fabrication  process. 

The  powder  is  pressed  (isostatically)  into  a  green  body  and  then  fired.  During 
cool-down,  it  is  possible  that  a  small  amount  (i.e.,  1  percent  or  less)  of  YSZ 
is  not  fully  stabilized  and  thus  around  2000°F,  the  structure  changes  from  cubic 
to  monoclinic  which  is  accompanied  by  a  7-10  percent  volume  change  of  the  mate¬ 
rial.  This  creates  high  local  stresses  which  remain  as  residual  stresses  when 
the  brick  eventually  cools  down  to  room  temperature.  Another  source  of  residual 
stresses  is  due  to  the  presence  of  glass  in  the  grain  boundaries  of  the  YSZ 
solid  solution.  Again,  during  cool  down,  the  viscous  nature  of  this  glass  will 
allow  stresses  to  redistribute  internally  and  eventually  to  result  in  locked  in 
residual  stresses  at  room  temperature.  These  residual  stresses  can  be  quite 
large  and  tend  to  accelerate  the  possibility  of  fracture  occurring  during  opera¬ 
tion  of  the  bricks  in  the  heater  bed. 

Locally  severe  thermal  stresses  would  be  associated  with  any  localized 
but  significant  perturbations  in  the  brick  matrix  temperature  field.  Such  local 
hot  spots  (or  cold  spots)  could  arise  from  a  variety  of  causes,  most  of  which 
would  be  associated  with  design  or  fabrication  imperfections.  For  example, 
local  increases  in  the  convective  heat  transfer  coefficient  due  to  hole  mis¬ 
alignments  at  the  brick  interfaces  (as  considered  as  part  of  Appendix  B)  would 
result  in  a  local  cold  spot  during  blowdown.  Similarly,  irregular J ties  at  the 
interface  between  the  outer  boundary  of  the  brick  matrix  and  the  heater  vessel 
interior  insulation  could  give  rise  to  stress  affecting  local  temperature  varia¬ 
tions. 
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Partial  destabilization  of  YSZ  could  possibly  occur  when  it  passes 
through  the  temperature  range  of  1900-2000*F.  This  causes  an  increase  in  the 
monoclinic  structure  (as  compared  to  the  stabilized  cubic  structure)  which  has 
an  accompanying  volume  change  (c.f.,  Ref.  3).  As  described  above,  this  local 
volume  change  gives  rise  to  quite  large  internal  stresses  which,  like  residual 
stresses,  can  be  quite  large. 

Micro-cracks  can  exist  in  the  bricks  due  to  the  fabrication  process.  As 
explained  above,  the  mixture  of  yttria  and  zirconia  may  not  be  completely  homo¬ 
geneous.  Consequently,  during  fabrication,  local  concentrations  of  unstabilized 
zirconia  will  experience  a  phase  change,  accompanied  by  the  volume  change,  and 
therefore  create  local  stresses  which  may  be  high  enough  to  create  local  cracks 
in  the  material.  This  phenomenon  is  slightly  different  than  the  residual  stress 
phenomenon  since  it  occurs  on  a  microscopic  scale  rather  than  a  macroscopic 
scale  (residual  stresses).  In  addition,  other  phenomena  also  cause  micro-cracks 
to  occur  in  the  fabricated  bricks.  The  result  of  these  micro-cracks  is  that 
they  act  as  fracture  triggering  mechanisms,  allowing  macro-cracks  to  initiate 
and  run  through  the  material,  thereby  producing  a  fracture  condition. 

All  of  the  above  listed  phenomena,  plus  others  not  listed  here,  have  the 
general  tendency  to  aggravate  the  stress  conditions  which  will  eventually  lead 
to  fracture  of  the  bricks.  Assuming  for  the  moment  that  all  of  these  other 
sources  give  rise  to  a  tensile  pre-stress  of  2000  psi,  then  for  the  example  of 
combined  stresses  shown  in  Figure  22,  the  ultimate  stress  needed  for  fracture 
would  only  be  6384  psi  (as  compared  to  8384  psi)  and  the  interaction  curve 
would  be  lowered  as  shown  by  the  long  dashed  line  in  Figure  22.  Now,  an  axial 
AT  *  150°F  and  a  radial  AT  =  40°F  would  suffice  to  cause  fracture,  as  indicated 
by  the  short  dashed  line  in  the  figure.  This  example  simply  illustrates  that 
these  other  stress  affecting  phenomena  tend  to  give  rise  to  stresses  which  can 
be  quite  large.  When  they  are  considered  in  conjunction  with  the  sources 
analytically  treated  in  this  report,  it  is  easily  demonstrated  that  certain  com¬ 
binations  of  all  these  sources  can  lead  to  fracture  conditions  within  the  bricks. 
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SECTION  4 


CONCLUSIONS 

1.  Relatively  precise  calculation  of  web  stresses  was  performed.  For  the 
conditions  considered,  the  maximum  predicted  tensile  stress  was  about  4000  psi. 
Recent  measurements  indicate  that  the  ultimate  strength  of  yttria  stabilized 
zirconia  is  about  8000  psi  at  this  condition.  It  is  concluded  that,  for  these 
conditions,  brick  fracture  is  highly  unlikely  due  to  web  stresses  alone,  parti¬ 
cularly  when  the  effect  of  the  heater  hydrostatic  pressure  is  considered. 

2.  Brick  normal  stresses  due  to  temperature  distribution  nonlinearities  in 
the  radial  and  axial  directions  were  computed.  The  results  were  correlated  in 
terms  of  the  relation 


°i 


-  6i 


E  a  AT 

TT““T 


r,0,Z 


where  the  6^  factors  were  defined  for  a  variety  of  conditions.  The  actual  mag¬ 
nitudes  of  temperature  nonlinearities  across  the  bricks  are  unknown,  so  these 
results  can  only  be  interpreted  in  a  parametric  fashion.  However,  an  approxi¬ 
mate  thermal  analysis  (Appendix  B)  indicates  that  radial  temperature  nonlineari¬ 
ties  of  the  order  100°R  are  possible  for  bricks  with  misaligned  holes.  This 
results  in  maximum  tensile  stresses  more  than  half  the  ultimate  strength. 


3.  It  is  unlikely  that  point  contact  loading  of  the  bricks  will  contribute 

significantly  to  brick  stresses. 


4.  The  principal  of  superposition  was  employed  to  assess  the  net  effect  on 
brick  stresses  of  combined  thermal  and  mechanical  loading.  Again,  the  magnitudes 
of  brick  net  temperature  distribution  nonlinearities  are  the  major  unknown.  An 
approximate  interaction  curve  was  constructed  which  considers  combined  loading 
and  gives  the  fracture  conditions  in  terms  of  required  radial  and  axial  tempera¬ 
ture  nonlinearities.  Limitations  on  the  use  of  this  curve  are  discussed  in 
Appe—  x  C. 


5.  The  brick  stresses  and  fracture  conditions  will  also  be  affected  by 
phenomena  which  are  not  amenable  to  theoretical  quantitative  analysis.  These 
include  residual  stresses,  incomplete  stabilization  or  destabilization,  micro¬ 
cracks,  and  local  hot  spots. 
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Uncertainties  relative  to  the  magnitudes  of  brick  temperature  distribu¬ 
tion  nonlinearities  and  the  factors  listed  in  Item  5  above  preclude  quantitative 
specification  of  heater  operation  modifications  or  erred  brick  design  revisions 
based  on  the  results  of  this  study. 


63 


■  as  ( Sa hi  i  r  PtrtiM  Wi 


'SraiTA-  T,%.' 


REFERENCES 


Hendrick,  W.  S.,  et  al.,  "Storage  Heater  Design  Study  for  the  Hypersonic 
True  Temperature  Tunnel,"  AEDC  TDR-64-48,  July  1964. 

Harney,  Major  D.  J.,  "Performance  Analysis  of  Large  Ceramic  Storage  Heaters," 
AEDC  TR-65-224 ,  December  1965. 

Hagford,  D.  E.,  DeCoursin,  D.  G.,  "Research  on  Storage  Heaters  for  High 
Temperature  Wind  Tunnels  -  Final  Report,"  AEDC  TR-71-258,  February  1972. 

Moyer,  C.  B.,  Davis,  J.  N.,  "Program  for  Analysis  of  Null  Point  Calorimeter 
Data,  Version  2  (PANDA  2)  Description  and  Users  Manual,"  Aerotherm  Report 
UM-72-29 ,  April  1972. 

Kreith,  F.,  Principles  of  Heat  Transfer,  Second  Ed.,  Int.  Texbook  Co.,  July 
1968,  Chapter  S. 

Schneider,  P.  J.,  Temperature  Response  Charts,  John  Wiley  and  Sons,  Chart 
36(c),  1963. 

Jacob,  M.,  Heat  Transfer,  Vol .  II,  John  Wiley  and  Sons,  1957,  Chapter  35. 

Moyer,  C.  B.,  " Ax i -Symmetric  Transient  Heating  and  Material  Ablation  Pro¬ 
gram  (ASTHMA)  Description  and  User's  Manual,"  Aerotherm  Corporation, 

Mountain  View,  California,  Aerotherm  Report  No.  68-27,  January,  1968. 

Engelhe,  w. ,  rersonal  Communication,  Southern  Research  Institute,  Birming¬ 
ham,  Alabama,  July  1971. 

Southern  Research  Institute,  "The  Thermal  Properties  of  an  Yttria  Stabilized 
Zirconia,"  Birmingham,  Alabama,  January  1972. 

Fahrenbacher ,  Major  L.  R. ,  personal  Communications,  ARL,  Wright  Patterson 
Air  Force  Base,  Ohio,  July  1971. 

Southern  Research  Institute,  "Preliminary  Mechanical  Property  Data,  Yttria 
Stabilized  Zirconica,"  Birmingham,  Alabama,  October  1971. 

Weiler,  F.  C.,  "DOASIS  Code  User's  Manual,"  Weiler  Research, Inc.  Report 
No.  1,  1971. 

Zienkiewicz ,  0.  C.,  The  Finite  Element  Method  in  Engineering  Science, 

McGraw  Hill,  London,  1971. 

Przemieniecki ,  J.  S.,  Theory  of  Matrix  Structural  Analysis,  McGraw  Hill,  . 

New  York,  1968.  \(>  ■ 


17.  Fung,  Y.  C. ,  Foundations  of  Solid  Mechanics,  Prentice-Hall,  Inc.  1965, 
Chapter  14 . 

18.  Eabel,  H.,  "EiaxiaJ  Fracture  Strength  of  Brittle  Materials,"  AFML-TR-66-51 , 
Air  Force  Materials  Laboratory,  Wright  Patterson  Air  Force  Base,  Ohio, 

May  1966. 


64 


LIST  OF  SYMBOLS 


A^  nodal  area,  node  i 

constant  in  Eq.  (A -6),  defined  by  Eg.  (A-7) 

constant  in  Eq.  (A-6) ,  defined  by  Eq.  (A-8) 

* 

b^  value  of  modified  by  recursion  relation  (A-15) 

c ^  constant  in  Eq.  (A-6),  defined  by  Eq.  (A-9) 

C  specific  heat 

P 

d^  constant  in  Eq.  (A-6),  defined  by  Eq.  (A-10) 

* 

d^  value  of  d^  modified  by  recursion  relation  (A-15) 

k  thermal  conductivity 

nodal  thickness,  node  i 
NL  denotes  last  node 

qcon<3  heat  conducted  into  body  away  from  heated  surface 

r  radial  distance  from  the  center  of  the  cylinder 

T  temperature 

Ta  initial  temperature 

Tq,Tw  surface  temperature 

u  thermal  conductance  defined  by  Eqs.  (A-ll)  ,  (A-12) , 

(A-13) 

a  constant  in  Eq.  (A-18)  defined  by  Eq.  (A-19) 

p,  constant  in  Eq.  (A-18)  defined  by  Eq.  (A-20) 

AO  finite  difference  time  st«.p 

ft  ^  i  **%«"* 

p  density 

t  nodal  time  constant,  Eqs.  (A-23) ,  (A— 2 4 ) 
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LIST  OP  SYMBOLS  (Concluded) 


Subscripts 

i 


denotes  node  i 


Superscripts 

n  denotes  at  time  step  n 
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APPENDIX  A 


IMPLICIT  ONE-DIMENSIONAL  FINITE  DIFFERENCE  CONDUCTION  SOLUTION 

The  solution  of  the  transient  heat  conduction  problem  for  a  variable- 
property,  one-dimensional  circular  cylinder  is  described. 

The  solution  uses  a  simple  implicit  finite  difference  approach.  The  dif¬ 
ferential  equation  for  T(rfe)  is 


—  ( kA  cC  A  il 

;.r  V  Or  )  P  SC 


( A- 1 ) 


with  the  boundary  conditions 


VT| 

r 

r=r 


=  0 


(A- 2) 


T I  =  T  (0)  =  T  (c)  (A- 3) 

! r=ri  ' 

where  TQ(f)  is  to  bo  determine^  from  the  system  energy  balance.  The  initial 
condition  is 


T(r,o)  =  T 

d 


(A- 4) 


where  T  is  a  constant. 

3 

In  the  difference  approach,  the  cylinder  is  broken  into  a  number  of  zones 
termed  nodes,  as  shown  in  Figure  A-l  . 

\  ( 

l 

!  •  • 

/ 

i= 1  i  =  2  . . .  i  «  NL 

Figure  A-l  Sketch  of  Nodal  C-r;d 
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The  temperature  is  defined  at  the  center  of  each  node  and  also  at  the  heated 
surface.  The  difference  analog  of  the  differential  equation  (A-l)  which  results 
from  taking  a  lumped  parameter  approach  is,  for  each  node  except  the  first  and 
last : 


(Tfl  -  Tt 


(»Vsi)  «i  - 


n+1  _  Tn+1 
li-l  i 

*i-i  4  *i 

'i-ih-i 


Tn+1  _  „n+l 
Ti  Ti+1 


(A-5) 


ci+lAi+l 


Here  T?  is  the  known  fold) temperature  of  node  i  at  time  0  (time  step  n) ,  and  T^ 
is  the  unknown  (new)  temperature  of  node  i  at  time  6  +  A6  (time  step  n  +  1) .  For 
the  last  (i  *  NL)  node,  the  last  term  in  the  right  hand  side  is  missing,  since  no 
heat  flux  leaves  the  insulated  back  wall.  For  the  first  node,  the 
quantity  is  missing. 

The  solution  procedure  requires  finding  all  the  "new"  nodal  temperatures 
T?+^  at  each  time  step.  A  new  time  step  ic  selected,  and  the  process  repeated  until 
the  final  time  is  reached.  For  solution  purposes.  Equation  (A-5)  is  written  as 


a.T?+}  +  b .  T*?+1  +  C.T?+} 
li-l  li  ii+l 


(A-6) 


where 


ai  =  ui 


(A-7) 


bL  =  -  ui  -  ui+1  -  pCp^i^A6  i  =  2,  ...  NL  (A— 8 ) 
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(A-12) 
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UNL+1  *  0  (A-13> 

The  entire  set  cf  nodal  equations  (A-6)  may  be  arranged  as  (dropping  the  super¬ 
script  for  temporary  convenience) . 


a, T  +  b,T,  + 
1  c  11 
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(A-14) 


a3T2  +  b3T3  +  C3T4 

V 


aNI*l‘NL-2  +  bNL-lTNL-l 


cnl-itml 


d 


3 


dKL-l 


aNLTn.-l 


+  b 


t:i.‘NL 


=  d 


NL 


Since  the  last  (lowest)  equation  has  only  two  unknowns,  it  may  be 
ir.ate  one  unknown  in  the  next  to  last  equation,  which  in  turn  may 
Thus  all  the  c,  terms  may  be  wiped  out  of  equations  (A-14)  in  the 
the  familiar  C-auss  reduction;  this  results  in  modified  b^  and 
to  the  follovinq  recur,  ion  relations; 
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After  this  pass,  the  first  equation  then  is 


Since  the  TACH  solution  procedure  is  implicit,  there  are  no  stability 
limits  on  the  choice  of  time  step.  However,  many  users  still  prefer  to  scale 
the  time  step  used  on  the  natural  stability  limit  of  the  nodal  system.  Therefore 
the  program  computes  the  natural  time  constant  for  all  nodes  at  each  time  step 


j  =  node  index 


(A-23) 


and  then  seeks  the  maximum  time  constant 


Tmax  =  max  l:t] 


(A-24 ) 


The  time  step  selected  is  then  t  if  t  does  not  exceed  the  user  input 

max  max 

maximum  time  step  or  fall  below  the  input  minimum  time  step. 

The  time  step  can  of  course  not  exceed  the  time  remaining  to  the  next 
output  time. 
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APPENDIX  B 

EFFECTS  OF  HOLE  MISALIGNMENTS 
IN  WIND  TUNNEL  HEATER  CORE  BRICKS 
ON  PRESSURE  DROP  AND  HEAT  TRANSFER 

by 

M.  D.  Jackson 
C.  A.  Powars 


LIST  OF  SYMBOLS 


A,  smallest  area  associated  with  sudden  expansion  or  sudden  contraction 

1  (ft2) 

A~  largest  area  associated  with  sudden  expansion  or  sudden  contraction 

2  (ft3) 

f  friction  factor 

G  mass  flux  (lbm/ft2sec) 

h  heat  transfer  coefficient  (Btu/hr-ft2-°F) 

K  nondimensional  head  loss  coefficient 

m  mass  flow  per  hole  (lbm/sec  hole) 

r  radius  (ft) 

R  gas  constant 

Subscripts 
avg  average 

cont  contraction 

exp  expansion 

fd  fully  developed 

max  maximum 

x  axial  distance 
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SECTION  B.l 


INTRODUCTION 

This  appendix  summarizes  the  results  of  a  brief  study  directed  toward 
estimating  the  effects  on  pressure  drop  and  heat  transfer  of  misalignments 
in  the  holes  through  ceramic  bricks  used  in  hypersonic  wind  tunnel  storage 
heaters. 

The  bricks  and  conditions  considered  correspond  to  the  A?TU  (Airbreathing 
Propulsion  Test  Unit)  facility  to  be  installed  at  AEDC.  These  conditions  are 
summarized  in  Table  B-l  below. 

TABLE  B-l 

WIND  TUNNEL  HEATER  CONDITIONS  CONSIDERED 

Nominal  pressure  *  3000  psi 
Nominal  flow  rate  =  1150  lbm/sec 

Initial  bed  temperature  distribution  is  sinusoidal  with 

Tbottom  "  1520°R  and  Ttop  *  4500°R 
Hole  diameter  (D)  =0.20  inches 

Hole  center  distance  to  diameter  ratio  (S/D)  =  1.47 
Heater  core  height  =  24  ft 

Yttria  Stabilized  Zirconia  bricks  in  top  12  ft  of  heater 
Alumina  bricks  in  bottom  12  ft  of  heater 
All  bricks  assu ned  to  be  6  inches  long. 

This  brief  study  was  motivated  by  the  fact  that  certain  Zirconia  bricks 
received  by  AEDC  for  use  in  a  pilot  heater  had  hole  misalignments  significantly 
greater  than  anticipated.  Quantitative  details  relative  to  the  hole  misalign¬ 
ments  are  available  from  Reference  B-l.  Three  cases  were  considered  here: 

1.  0%  misalignment  (i.e.,  no  misalignment;  the  baseline  case) 

* 

2.  14%  misalignment  (corresponds  to  an  approximate  average  misalignment) 

* 

3.  30%  misalignment  (corresponds  approximately  to  worst  case  misalignment 

considering  keying  tolerance  and  hole  misalignments) 

J - 

Percent  misalignments  are  defined  on  an  area  basis;  blocked  area/area  originally 
open  to  flow. 


The  effects  of  the  hole  misalignments  on  fluid  friction  and  pressure  drop 
are  discussed  in  Section  B.2.  Estimates  of  the  heater  core  pressures  differentials 
for  a  constant  flow  rate,  and  of  the  flow  rates  through  the  bricks  for  a  constant 
heater  core  pressure  drop,  are  given.  The  effects  of  the  hole  misalignments  on 
heat  transfer  from  the  heater  core  to  the  fluid  being  heated  are  discussed  in 
Section  B.3.  A  quantitative  calculation  of  the  brick  temperature  distribution 
histories  was  beyond  the  scope  of  this  study,  although  certain  limiting  cases  are 
considered.  Conclusions  are  drawn  and  recommendations  are  made  in  Section  B.4. 


SECTION  B.2 


EFFECTS  ON  PRESSURE  DROP  AND  FLOW  RATE 

The  effects  of  the  brick  hole  misalignments  on  the  heater  pressure  drop 
and  flow  rates  are  considered  in  the  following  subsections. 

B.2.1  GOVERNING  EQUATIONS 

The  pressure  drop  through  aligned  brick  holes  may  be  easily  calculated  from 
a  mass  and  momentum  balance.  For  a  constant  mass  flow  rate  and  a  variable  density 
fluid  the  momentum  equation  expresses  a  balance  between  acceleration,  pressure 
force,  and  wall  shear.  The  wall  shear  for  turbulent  flow  in  pipes  may  be  defined  in 
terms  of  a  friction  factor  and  substituted  into  the  momentum  equation.  Also,  us¬ 
ing  continuity  and  the  equation  of  state  for  an  ideal  gas,  one  can  express  the 
momentum  equation  in  terms  of  temperature,  mass  flux  and  pressure  as  follows: 

§1 +  tt(?)  +  g2r  3x(?)  "  0  (B“1) 

where 

G  =  mass  flux 
f  =  friction  factor 
r  =  radius 
R  *  gas  constant 

Equation  (B-l)  is  the  governing  equation,  valid  for  one  dimensional  steady  turbu¬ 
lent  flow  and  can  be  solved  for  dP/dx  if  the  temperature  distribution  is  given. 

Equation  (B-l)  may  be  simplified  if  the  temperature  gradient  is  much 
greater  than  the  pressure  gradient.  This  is  the  case  for  the  conditions  of 
interest  here;  it  can  be  shown  that 


to  within  1%  error. 


d  /t\  1  dT 

cHF  \^p  j  '  p  3x 

Hence  Equation  (B-l)  reduces  to 
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(B-2) 


Misalignment  of  holes  at  brick  interfaces  will  cause  the  flow  to  separate, 
Ideally  the  pressure  loss  can  be  calculated  from  a  mass  and  momentum  balance. 


TH 
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However,  for  flows  involving  separation  it  is  convenient  to  express  the  pressure 
loss  in  terms  of  a  nondimensional  head  loss  (K)  where  K  is  defined  as 


K 


-  _A!Vc 
1/2  V  2 


(B-3) 


For  turbulent  flows  in  which  separation  occurs  K  is  to  a  good  approximation 
dependent  only  upon  boundary  geometry  or  the  amount  of  separation  caused  by  a 
particular  boundary  geometry.  Thus  K  can  be  estimated  from  existing  literature 
values  for  systems  with  similar  geometry. 

Equation  (B-3)  is  valid  for  steady  incompressible  flow.  Since  the  region 
of  separation  for  the  mass  fluxes  encountered  in  this  problem  is  sufficiently 
small  compared  to  nominal  brick  length,  one  can  assume  that  Equation  (3-3)  is 
valid  at  the  brick  interface.  Solvirg  for  AP  and  using  the  ideal  gas  equation 
of  state,  Equation  (B-3)  can  be  expressed  as  follows: 


K  G2TR 
2  P 


(B-4) 


The  determination  of  K  will  be  considered  in  Section  B.2.2. 

Since  Equation  (B-4)  is  applicable  only  at  discrete  sections,  at  which 
T,  P  are  dependent  on  upstream  conditions,  it  is  necessary  to  rewrite  Equation 
( B—  2 )  as 
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( B—  5 ) 


The  total  pressure  drop  across  the  heater  bed  may  be  calculated  by  progressing 
downstream  through  a  hole  applying  Equation  (B-5)  along  the  length  of  a  brick  and 
Equation  (B-4)  at  the  brick  interface. 

B.2.2  ESTIMATION  OF  HEAD  LOSS  AT  BRICK  INTERFACE 

As  mentioned  in  Section  B.2.1,  K  is  a  convenient  correlation  which  can 

be  used  to  determine  the  pressure  drop  duo  to  misalignment .  Figure  B-l  shows  the 

functional  relationship  between  K  and  area  ratio  for  a  sudden  expansion  and  sudden 

contraction  (note  that  .or  0.5  *■  A./A-  <  1.0,  K  -  K  .  )  .  For  misaligned 

—  l  i  —  exp  cont 

holes  the  resistance  to  flow  is  area ter  than  a  sudden  expansion  and  less  than  a 
sudden  contraction.  Hence  if  an  equivalent  area  ratio  can  be  determined  at  brick 
interfaces  and  provided  this  area  ratio  is  between  0. 5-1.0  a  reasonable  estimate 
can  be  made  cf  K. 
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For  turbulent  flow  K  is  a  function  only  of  system  geometry,  hence  one 
can  define  the  area  ratio  in  terms  of  equivalent  flow  blockage  for  a  sudden 
contraction  or  conversely  in  terms  of  equivalent  flow  expansion  for  a  sudden 
expansion.  As  pointed  out  in  the  introduction  the  two  cases  of  area  misalign¬ 
ment  of  interest  were  14%  and  30%.  Using  the  above  formulation  for  area  ratio 
one  obtains  K  =  .026  and  K  -  .10  for  14%  (\/*2  m  ,86)  and  30%  (*]/A2  *  *70) 

misalignment  respectively.  These  values  of  the  nondimensional  head  loss  are 
the  average  of  sudden  expansion  and  sudden  contraction  at  the  given  area  ratio. 

The  head  loss  at  the  brick  interfaces  was  also  estimated  by  assuming 
that  K  is  proportional  to  the  volume  of  separation,  and  by  applying  a  momentum 
balance  while  assuming  all  energy  dissipation  occurs  in  the  region  of  expanding 
flow.  These  alternative  approaches  yielded  values  of  K  which  were  in  essential 
agreement  with  the  values  given  above. 

B.2.3  CALCULATION  OF  BED  PRESSURE  DROP  FOR  0%,  14%,  and  30%  HOLE  MISALIGNMENT 

WITH  CONSTANT  FLOW  RATE 

Section  B.2.1  briefly  summarized  the  equations  used  to  estimate  the  total 
bed  pressure  drop  through  the  heater  bricks.  Using  the  values  of  K  found  in 
Section  B.2.2  the  total  pressure  drop  can  be  computed  by  progressing  through  the 
heater  bed  using  Equation  (B-5)  along  a  brick  length  and  applying  Equation  (B-4)  at  i 
brick  interface. 

This  calculation  was  carried  out  for  a  total  heater  flow  rate  of  1150 
lbm/sec,  surface  roughness  of  100  u  inches  and  for  a  temperature  distribution 
corresponding  to  10  seconds  from  the  initial  time.  The  results  are  shown  in 
Tables  B-2  and  B-3  for  14%  and  30%  area  offset.  For  zero  percent  misalignment, 
the  pressure  drop  across  the  bed  is  4.8  psi  for  these  conditions.  Hole  mis¬ 
alignments  of  14%  and  30%  increased  the  total  pressure  drop  to  6.7  psi  and 
12.3  psi  respectively.  Note  that  at  30%  misalignment  the  total  pressure  drop 
is  still  only  0.4%  of  the  entrance  pressure. 

B.2.4  CALCULATION  OF  FLOW  RATES  FOR  0%,  14%,  AND  30%  HOLE  MISALIGNMENTS  WITH 

CONSTANT  PRESSURE  DROP 

In  Section  B.2.3  the  pressure  drop  across  the  heater  was  estimated  assuming 
a  constant  heater  flow  rate  (li’^O  lbm/sec)  for  hole  misalignments  of  0%,  14%, 
and  30%.  In  reality,  the  misalignment  of  the  holes  is  not  uniform  and  varies 
from  near  0%  at  the  brick  center  to  higher  values  closer  to  the  edge  of  the  brick. 
Since  the  pressure  drop  across  the  bed  must  be  the  same  for  all  holes,  misaligned 
holes  will  have  lower  flow  rates  than  aligned  holes.  This  variation  in  flow 
rates  is  estimated  in  this  section. 


80 


iIAL  LkJCATIPK 

FRICT*THfcRM  HRLSS  OR(JP  0FF 

PRESS  DR0P 

78TAL  PRESSURE 

,60 

,03812 

,00771 

2999,95419 

1  ,CG 

.04356 

,01018 

2999,90048 

1,60 

,04738 

.01193 

2999,04122 

2,(0 

, 0609C 

,01353 

2999,77662 

2,60 

.05473 

,01527 

2999,70686 

3,00 

,05856 

,01702 

2999,63129 

3 , 50 

,06272 

,01891 

2999,54960 

4,00 

.06691 

,02036 

2999,46344 

4,60 

,06943 

,02196 

2999,37207 

6,00 

,07295 

,02357 

2999,2755/ 

5.6U 

,07646 

,02617 

2999, 17J98 

6, CO 

,08142 

,02742 

2999,06519 

6,60 

,06493 

,02902 

2998,95123 

7,  CO 

,074/5 

,02924 

2990,84/26 

7,60 

,07827 

,03084 

2998,73819 

8,00 

,06146 

,03230 

2990,62445 

8,60 

.08466 

, 03376 

2996 , 506C4 

9.  CO 

,06722 

,03492 

2998,38391 

9,60 

.05042 

, 03638 

2998,25714 

10. CO 

.09362 

,03783 

2998,12573 

10,60 

,09618 

,03900 

2997,99060 

11, CO 

,09664 

,04017 

2997,85382 

1 1 ,5C 

,09984 

,04162 

2997,71240 

12,00 

,10208 

,04264 

2997,56/69 

12.60 

,10433 

,04367 

299/, 41974 

13,00 

.10669 

,04483 

2997,26804 

13,60 

,10913 

,04585 

2997,11310 

14,00 

,11137 

,04687 

2996,95490 

14,50 

.11362 

,0479U 

2996,79340 

15, CO 

,11666 

,04892 

2990,62866 

16,50 

. 11779 

,04979 

2996,46112 

16,00 

,12035 

,05096 

2996,28966 

16,50 

,12196 

,05169 

2996,11624 

17,  CO 

,12420 

,05272 

2995,93936 

17,50 

.12645 

,05374 

2995,75922 

18, CO 

,12596 

,05461 

2995,57067 

18,50 

,12789 

.05549 

2995,39632 

19,00 

.12962 

,05637 

2995,20917 

19.50 

. 13174 

,05/25 

2995,02023 

20, CO 

, 1 3335 

,0579b 

2994.82895 

20,60 

.  *.3560 

,05900 

2994,63437 

21,00 

. i361C 

,05944 

2994,43966 

21,50 

. 136J9 

,0',00J 

2994,24347 

22, Oo 

,13/68 

,  1)6062 

2964,04520 

22,50 

, 13897 

,06120 

2993.845C6 

23,00 

, 1 4058 

,06194 

2993,64258 

23,50 

.14155 

62  «8 

2993,43669 

24, CO 

. 1 4220 

,00000 

2993,29651 

TfcTAl  FKKT1PN* 

4,0210  HSi  TeT*i.  ph-sft 

*  1,8640 

PS1 

TABLE  B-2 

A.XIAL  PRESSURE  DROP  ACROSS  HEATER  BED,  14%  MISALIGNMENT 


81 


1*1  10CAT1PN 
,  bo 
l  .ro 
i  .ho 
2 .  C  0 
2, bo 
'3. CO 

3. 50 

*.ro 

4.50 

5. (u 
b.bo 
p.co 

6. M 

7. ro 
7.  bo 
*.C( 
e.bt 
o.rc 
o.bo 

10.  CO 
It'. bo 

11. ro 

it  ,bO 

12. ro 
12.  bo 
iJ.ro 

1  J,5u 

14. (0 
14.SII 

15.  CO 
Ib.bt 

IO,  CO- 

16.  50 

17. C0 

17.50 
lft.r-0 

IP.  bo 
10. to 
1^.50 
20.(0 

20.50 
21,(0 
21  .bo 
2?. r  t. 

22.50 
2.1, CO 
21, bo 
24.ro 

t»  i*i  fKii  im« 


(•PIC  !*TH£RH  P  R  L  5  S  0*66  *** 

.03*12 
.04  5 55 
,047j8 
.Otiogr 
.0547.3 
, 05857 
.0627? 

.06592 
,06644 
.07296 
."7647 
.Obi 4J 
,0H49b 
,07477 
."7*29 
, Ob  146 
,"H4o9 
."6725 
.09046 
,  09366 
,0962.3 
,09669 
,09990 
.10214 
.  t  04.39 
.10*96 
.10921 
.11146 
.11  371 
,11596 
.11790 
.12047 
,12211* 

.124J4 
.12656 
.12611 
.12*05 
.  1299* 

.1  J192 
,  1.3  354 
,1 J6*C 
.1 J5J1 
,  1.3661 
. 1 J/9C 
,1 J920 
.  1  4i’6.3 
,14161 
,  1  4247 

4,625.3  PM  1  v*  1  A l.  Pf3 


PPF.SS  0H14P 

T8T*L  PPtSOUHt 

.03063 

2999.93106 

.04072 

2999,84660 

,04/70 

2999,75174 

,05411 

2999,64679 

,06109 

2999,53101 

,06607 

2999,40439 

.07564 

2999,268(5 

,06146 

2999,11868 

,06787 

2998,96140 

,09427 

2996,79419 

,10066 

?096, 61/0  7 

.10971 

2998,42590 

.11612 

2998,22495 

,11700 

2998,03323 

.12J41 

299/, 63157 

,12924 

2997,62068 

.  1.3507 

2997,40115 

.13974 

290/,  17419 

,14557 

2996.93817 

,  1 5  i  4  1 

2996.09312 

,15608 

2990,44063 

,1607b 

2990. 1P341 

,16659 

2995.9169b 

,l706e 

29C5. 04417 

.17477 

2995.36305 

,17945 

2995,0706/ 

,18355 

2994,78595 

,16764 

2994,460*0 

.19174 

2994.18J43 

,19564 

2993.86960 

,1993b 

2993,5524 j 

,20405 

2993,22794 

,20698 

2992.89690 

,21109 

2992.56351 

,21519 

2992.22174 

.21872 

2991,87095 

,2222 4 

2991,520/0 

,22577 

2991,17099 

,229.30 

2990,60961 

.23224 

299L.444C  0 

,23635 

2990,07196 

,23814 

29P9, 69855 

,24050 

29»9, .*2147 

,2428/ 

7986,940/J 

,24523 

2968,55634 

.24816 

2968,16/30 

.24997 

29*7.77563 

,1)0000 

29*7,63318 

T »  7.5430 

PSI 
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TABLE  B-3 


AXIAL  PRESSURE  DROP  ACROSS  HEATER  BED, 


30%  MISALIGNMENT 
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To  obtain  a  bimple  relation  for  calculating  the  mass  flow  rate  per  hole 
at  a  given  pressure  drop  and  degree  of  misalignment,  it  was  assumed  that  the  total 
pressure  droo  could  be  expressed  in  the  functional  form  of  Equation  (B-4) 

This  amounts  to  defining  an  effective  net  head  loss  for  0%,  14%,  and  30% 
misalignment.  Using  Equation  (B— 3)  the  following  ratio  may  be  formed: 


mo%,  14%,  30%  =  (*0144)\a: 


across  bed 
30%,  14%,  30% 


(B-6) 


where  m  is  the  mass  flow  per  hole  and  has  units  of  lbm/sec.  The  fi P's  refer  to 
those  calculated  in  Section  B.2.3  for  total  mass  flow  rate  of  1150  lbm/sec. 

From  Equation  (B-6),  the  following  mass  flows  rates  per  hole  were  computed: 

•  AP=. ,  =  4.8  psi 

across  bed  r 

m  (0%)  =  .0144  ibm/sec-hole 
m  (14%)  =  .0122  lbm/sec-hole 

m  (30%)  =  .0086  lbm/sec-hole 

•  -iP  .  .  -  6.7  psi 

across  bed  r 

m  (0%)  =  .0170  lbm/sec-hole 
m  (14%)  =  .0144  lbm/sec-hole 
m  (30%  =  .0104  lbm/sec-hole 

•  l  _  j  =  12.8  psi 
across  bed 

m  (0%)  =■  .0231  lbm/sec-hole 
m  (14%)  =  ,0199  lbm/sec-hole 

m  (30%)  =  .0144  lbm/sec-hole 

The  above  calculated  mass  flow  rates  will  be  used  in  Section  B.3.2  in 
estimating  the  magnitude  of  the  temperature  gradient  between  holes. 
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SECTION  B. 3 


EFFECTS  ON  HEAT  TRANSFER 

One  of  the  primary  reasons  for  accessing  the  thermal  response  of  the  heater 
core  bricks  is  to  predict  the  severity  of  brick  thermal  stresses.  Hole  misalign¬ 
ments  in  heater  cr.e  bricks  could  augment  thermal  stresses  via  two  mechanisms: 

1.  hole  misalignment  causes  the  flow  to  separate  and  reattach  resulting 
in  a  substantial  axial  variation  in  the  local  heat  transfer  rates;  anc 

2.  hole  misalignment  increases  the  flow  resistance,  reducing  the  mass 
flow  rate  and  hence  the  heat  transfer  such  that  significant  radial 
temperature  gradients  might  exist  across  bricks  having  variable 
hole  misalignments. 

The  objective  of  this  section  is  to  evaluate  the  importance  of  the  above 
mechanisms.  Section  B.3.1  discusses  the  effects  cr  separation  and  Section  B.3.2 
will  deal  with  temperature  gradients  between  holes.  The  heat  transfer  physics 
associated  with  the  above  mechanisms  are  of  sufficient  complexity  that  only 
approximate  estimates  are  made. 

B.3.1  HEAT  TRANSFER  IN  SEPARATED  FLOWS 

Separated  flows  cause  a  substantial  axial  variation  in  the  local  transfer 
coefficient,  nx,  and  hence  in  the  local  heating  rates.  Large  axial  variations 
in  heating  rates  could  cause  material  failure  due  to  thermal  stressing.  The 
fluid  mechanics  of  separated  flows  are  highly  complex  and  reasonable  analytical 
methods  for  predicting  the  flow  field  do  not  exist.  Therefore  it  is  not  possible 
to  analytically  predict  hx#  and  one  must  rely  on  heat  transfer  data. 

Typically  in  separated  flows  the  axial  variation  of  hx  follows  a  bell 
shaped  curve.  That  is  at  the  point  of  separation  hx  is  at  a  minimum  below  that 
which  exists  at  fully  developed  conditions  and  rises  as  one  proceeds  downstream 
toward  the  reattachment  point.  At  the  point  of  reattachment  hx  reaches  a 
maximum  and  then  decays  to  a  value  consistent  with  fully  developed  flow. 

The  data  of  Zemanick  and  Dougall  (Reference  B-2)  on  the  heat  transfer 
rates  downstream  of  an  abrupt  circular  channel  expansion  represents  tne  best 
data  available  as  far  as  being  geometrically  similar  to  the  problem  of  misaligned 
holes.  From  their  data  the  ratio  of  the  maximum  local  heat  transfer  coefficient 
to  the  local  fully  developed  heat  transfer  coefficient  can  be  estimated. 
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•  14%  Area  Blockage,  hmax/hfd  = 

•  30%  Area  Blockage,  h /h,,  ~  3.2 

max  fa 

The  maximum  hx  occurs  at  the  point  of  reattachment,  which  is  dependent  on  system 
geometry.  Zemanick  and  Dougall  data  indicate  that  for  A^/Aj  *  .70  the  point  of 
reattachment  is  approximately  0.1  inches  downstream  of  the  sudden  expansion. 

Hence  hx  changes  by  a  factor  of  3  in  0.1  inches.  This  large  variation  could  give 
rise  to  uery  high  localized  temperature  gradients  and  thermal  stresses. 

Separation  also  has  the  effect  of  increasing  the  overall  average  heat 
transfer  coefficient.  While  this  overall  increase  in  h  is  not  calculated  here, 
three  approaches  could  be  taken  to  estimate  this  increase: 

•  Integrating  hx  over  the  length  of  heater  bed,  where  hx  could  be 
estimated  from  the  data  of  Zemanick  and  Dougall. 

•  Estimating  an  equivalent  friction  factor  due  to  misalignment 
(analysis  similar  to  Section  B.2  4)  and  using  correlations  which  have 
been  developed  for  heat  transfer  on  roughened  surfaces  (see  for  ex¬ 
ample,  Reference  B-3). 

•  Utilizing  a  simple  bench-scale  apparatus  (e.g.,  misaligned  copper 
tubes  fitted  with  thermocouples)  to  measure  h  experimentally 

B.3.2  EFFECTS  OF  HOLE  MISALIGNMENTS  ON  HEAT  TRANSFER 

As  pointed  out  in  Section  B.2. 4,  a  decrease  in  the  mass  flow  rate  due  to 
misalignment  could  cause  a  substantial  radial  temperature  gradient  between  holes 
having  different  mass  flows.  This  problem  is  difficult  to  handle  analytically 
due  to  the  heat  transfer  coupling  between  holes.  If  one  assumes  that  the  holes 
are  independent,  i.e.,  no  cross  conduction,  and  tnat  misalignment  does  not 
increase  tne  overall  heat  transfer  coefficient,  then  an  upper  bound  on  the  tempera¬ 
ture  gradient  between  holes  can  be  computed. 

The  TACH  code  (described  in  the  main  text)  was  utilized  to  predict  the 
thermal  response  of  holes  with  different  mass  flow  rates  due  to  misalignments. 

In  Section  B.2. 4  mass  flow  rates  were  calculated  for  0%,  14%,  and  30% 
hole  misalignments  at  a  given  pressure  drop  across  heater  bed.  These  values  of 
mass  flow  rates  were  then  used  in  the  thermal  response  code.  Figures  B-2  and 
B-3  illustrate  the  effect  of  different  mass  flow  rates  on  the  heater  thermal 
response,  subject  to  the  assumptions  discussed.  As  can  be  seen  from  these 
figures,  the  temperature  difference  initially  increase  with  time. 
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The  maximum  temperature  difference  encountered  at  shut  down  time  is  of 
primary  importance.  Since  the  overall  heat  transfer  coefficient  for  turbulent 
flow  in  pipes  is  proportional  to  the  Reynolds  number  raised  to  the  .  '  power, 
the  largest  temperature  differences  will  occur  at  the  lower  mass  flow  rates  or 
at  the  lowest  pressure  drop  across  heater  bed.  Figure  B-4  illustrates  the  exit 
temperature  difference  for  AP  =  4.8  psi.  As  can  be  seen,  two  cases  are  consid¬ 
ered:  A1  based  on  0%  and  14%  misalignment  and  AT  based  on  0%  and  30%  misalign¬ 
ment.  Shut-down  time  was  determined  when  the  exit  temperature  for  0%  misalignment 
dropped  by  10%.  For  convenience  the  maximum  temperature  differences  at  shut 
down  time  are  summarized  below: 

•  AT  (14%  misalignment  -  0%  misalignment)  =  60°R 

•  AT  (30%  misalignment  -  0%  misalignment)  =  290°R 

Using  the  above  values  for  the  maximum  temperature  differences  at  shut 
down,  one  can  make  a  simple  calculation  to  determine  the  magnitude  of  conduction 
across  the  holes.  It  is  found  from  this  calculation  that  conduction  across 
holes  is  the  same  order  of  magnitude  as  convection  across  brick-air  interface. 
Therefore  it  can  be  concluded  that  the  estimated  maximum  temperature  differences 
are  overestimates,  i.e.,  such  a  large  temperature  difference  would  not  occur 
since  the  conduction  heat  transfer  coupling  would  tend  to  dissipate  this 
temperature  difference.  An  analytical  computation  of  this  multiply  coupled 
heat  transfer  problem  would  be  relatively  complex  and  was  beyond  the  scope 
of  the  preliminary  study. 
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SECTION  B.4 

SUMMARY  AND  CONCLUSIONS 

Pressure  Drop  and  Flow  Rate 

•  For  an  average  misalignment  of  14%  the  pressure  drop  across  the  heater 
bed  will  increase  from  4.8  psi  to  about  6.7  psi. 

•  Pressure  drop  will  be  less  than  1%  of  entrance  pressure  (3000  psi) 
for  all  cases  of  misalignment. 

•  For  heater  bed  pressure  drops  of  4.8  psi,  14%  and  30%  misalignments 
will  reduce  the  mass  flow  rates  per  hole  by  15%  and  38%  respectively. 

Heat  Transfer 

•  Due  to  flow  separation  and  reattachment,  the  local  heat  transfer 
coefficient  (hx)  will  vary  by  a  factor  of  about  2  or  3  within  0.1 
inches  of  the  brick  interfaces  for  misaligned  holes. 

•  The  influence  of  this  variation  in  hx  on  thermal  stresses  should  be 
estimated  analytically.  If  the  thermal  stress  calculation  shows 
that  a  variation  in  hx  is  critical,  the  hx  variation  could  be 
defined  experimentally  using  a  simple  bench  scale  apparatus. 

•  Heat  transfer  analyses  indicated  that  an  upper  limit  of  the  temperature 
gradient  across  a  brick  with  variable  hole  misalignments  is  about 

60°R  for  0  to  14%  misalignment  and  300°R  for  0  to  30%  misalignment. 
These  calculations  neglected  cross-conduction  and  heat  transfer 
coefficient  increase  for  misaligned  holes.  In  view  of  these  conserva¬ 
tive  assumptions,  it  is  estimated  that  the  actual  temperature 
gradients  across  the  bricks  will  be  less  than  half  this  large. 
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The  equivalent  cylinder  model  used  to  evaluate  the  hypothetical  non¬ 
linear  temperature  gradients  within  the  bricks  was  only  an  approximation  to  the 
real  brick's  geometry.  Since  the  real  nonlinear  gradients  were  not  known  ex¬ 
plicitly,  this  2-dimensional  model  was  used  for  economic  reasons,  however,  as 
shown  by  the  in-plane  analyses,  the  hexagonal  geometry  does  produce  stresses 
which  deviate  from  equivalent  cylinder  geometries  as  demonstrated  by  the  shear 

stress  r  (cf.  text  Figures  13e  and  14d).  In  a  cylindrical  geometry,  these 
xy 

shear  stresses  would  be  zero.  Consequently,  a  more  accurate  3-dimensional 
model  of  the  hexagonal  brick,  like  the  one  shown  in  Figure  C-l  would  yield  more 
accurate  results  for  the  nonlinear  temperature  gradient  loading  cases.  Since 
the  stresses  from  this  source  of  loading  proved  significant  in  contributing 
towards  the  total  stresses  which  could  produce  fracture,  a  more  accuraue  deter¬ 
mination  of  the  stress  state  within  a  hexagonal  brick  3-dimensior.al  model  would 
be  useful.  However,  the  cost  of  running  any  3-dimensional  model  is  high  and 
this  analysis  should  not  be  undertaken  without  knowing  explicitly  the  aegree  of 
nonlinear  temperature  gradients  which  could  exist  within  the  bricks.  Once  an 
accurate  determination  of  these  nonlinear  gradients  is  made,  then  a  3-dimensional 
stress  analysis  utilizing  a  model  similar  to  the  one  shown  in  Figure  C-l  would 
yield  results  which  are  more  accurate.  These  results  could  then  be  used  in  con¬ 
junction  with  the  in-plane  analyses  (which  are  quite  accurate)  to  produce  a 
better  interaction  curve  like  the  one  shown  in  Figure  22.  Since  the  interac¬ 
tion  curve  shown  in  Figure  22  was  based  upon  the  2-dimensional  equivalent 
cylinder  model,  it  should  only  be  used  in  a  qualitative  manner. 

Another  reason  for  only  using  this  curve  in  a  qualitative  manner  is  that 
the  fracture  criterion  used  did  not  take  into  account  the  shear  stresses  which 
were  present  within  the  bricks.  When  examining  the  in-plane  solutions,  one 
notices  that  fairly  significant  shear  stresses  can  develop  (cf.  text  Figures 
13e  and  14d).  In  fact,  these  shear  stresses  were  approximately  40  percent  of 
the  maximum  tensile  stresses  within  the  brick.  When  postulating  a  fracture 
criterion  for  yttria  stabilized  zirconia,  one  immediately  notices  that  the  de- 
formational  behavior  of  YSZ  and  other  ceramics  falls  into  the  class  of  materials 
which  are  elastic-brittle,  similar  to  concrete  and  rocks.  The  literature  abounds 
with  different  failure  criteria  for  this  class  of  materials  but  most  of  it 
centers  around  the  well  known  Coulomb-Mohr  failure  cri*-  rjon  in  one  form  or 
another.  A  typical  failure  surface  of  this  type  is  shown  in  Figure  C-2 ,  where 
Mohr's  theory,  represented  by  the  solid  dark  line,  and  Coulomb's  theory,  re¬ 
presented  by  the  dashed  line,  are  both  illustrated.  The  Mohr  theory  assumes 
that  the  limiting  shear  stress  is  t  =  f(o),  i.e.,  some  function  of  the  nox.i.al 
stress  acting  upon  a  small  element  of  material.  One  such  possible  functional 
relationship  would  be 
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t  2  =  4o  (o.  -  a  ) 
t  t  n 


where  =  ur.iaxial  tensile  strength,  =  normal  compressive  stress  and  t  = 
shear  rupture  stress. 

In  contrast,  a  typical  Coulomb  criterion  would  be 

T  =  T  +  U  ff 

o  n 

where  tq  =  intrinsic  shear  resistance  of  the  material  and  u  =  coefficient  of 
internal  friction,  and  t  and  on  are  the  same  as  above. 

Without  going  into  more  detail  (cf.  Reference  C-l  for  further  details), 
one  can  notice  that  the  failure  of  elastic-brittle  materials  can  depend  upon  more 
than  just  the  uniaxial  tensile  strength  of  the  material.  In  fact,  even  under 
compressive  loading  conditions,  when  high  shear  stresses  are  present,  failure  will 
occur  as  shown  by  the  Coulomb-Mohr  failure  envelopes  in  Figure  C-2.  Such  a  con¬ 
dition  exists  as  shown  for  the  point  contact  loading  analysis  as  shown  in  text 
Figures  20c  and  20e.  Although  these  stresses  were  not  severe  enough  to  cause 
failure,  the  type  of  stresses  present  were  compresive  normal  stresses  accompanied 
by  shear  stresses. 

The  main  point  to  be  made  here  is  the  following.  In  the  normal  stress 
"tensile"  region  (cf.  Figure  C-2),  when  shear  stresses  are  present,  the  magnitude 
of  the  tensile  stresses  needed  to  cause  fracture  are  less  than  the  uniaxial  ten¬ 
sile  strength  of  the  material.  When  the  shear  stresses  are  zero,  then  the  ten¬ 
sile  stresses  must  reach  the  uniaxial  tensile  strength  in  order  to  cause  fracture. 
Since  the  shear  stresses  were  neglected  in  the  combined  stress  analysis,  the 
interaction  curve  shown  in  text  Figure  22  represents  an  unconservative  estimate 
of  conditions  necessary  to  cause  fracture,  i.e.,  fracture  could  occur  for  values 
of  combined  stresses  less  than  that  indicated  by  this  interaction  curve.  There¬ 
fore,  as  stated  before  this  interaction  curve  should  only  be  used  in  a  qualitative 
manner,  and  not  as  a  design  curve. 
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